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ABSTRACT 

A  mixed-integer  linear  programming  formulation  is  developed  for 
minimizing  delay  to  traffic  in  a  signal  controlled  road  network.  Off¬ 
sets,  splits  of  green  time  and  a  conuion  cycle  time  for  the  network  are 
considered  as  decision  variables  simultaneously.  The  traffic  flow  pat¬ 
tern  is  modeled  as  a  periodic  platoon,  and  a  link  performance  function 
is  derived  in  the  form  of  a  piecewise  linear  convex  surface  representing 
the  delay  incurred  by  these  platoons.  Stochastic  effects  are  accounted 
for  by  a  saturation  deterrence  function  representing  the  expected  over¬ 
flow  queue  on  each  link  and  are  included  as  an  additive  component  in  the 
objective  function.  Computational  results,  using  the  MPSX  system,  are 
given  for  an  arterial  with  11  signals  in  Waltham,  Massachusetts,  and  a 
portion  of  the  UTCS  network  in  Washington,  D.C.  containing  20  nodes,  63 
links  and  21  loops. 
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The  last  fifteen  years  have  seen  a  gradual  infusion  of  computer 
technology  and  operations  research  methods  into  traffic  signal 
systems.  This  has  manifested  itself  in  three  main  ways: 

•*  in  offline  methods  for  calculating  signal  settings 

•  in  "solid-state"  digital  control  hardware,  and 

•  in  real-time  closed-loop  control  systems  using  both  on-line 
traffic  detection  and  calculation  of  settings. 

We  shall  describe  each  briefly  so  as  to  relate  our  work  to  the 
field  in  general. 

Offline  analytic  methods  may  be  considered  to  start  with  the 
works  of  Cl ay ton [5]  and  Webster  [40]  in  Britain,  for  determining 

qreen  SDlits  and  cycle  time  at  an  isolated  intersection.  For  coordinated 
signals  modern  research  begins  with  Morgan  and  Little's  rigorous 
maximization  of  bandwidth  on  an  artery[28]*  Little[24]  subsequently 
generalized  and  extended  the  methods  to  networks.  SIGOP  [37] 
introduced  a  quadratic  link  performance  function  and  a  new,  but  not 
global,  method  of  optimization.  Hillier  [15]  and  Allsop  [1] 
developed  the  Combination  Method.  This  permits  the  use  of  an 
arbitrary  link  performance  function  and  produces  global  optimal 
offsets  for  fixed  values  of  green  splits  and  cycle  time.  Gartner[8] 
formulated  the  problem  in  terms  of  dynamic  programing  optimization 
and  developed  an  efficient  network  partitioning  algorithm.  A 
different  direction  was  taken  by  Robertson  in  TRANSYT  [34]. 

Whereas  the  link  performance  function  approach  assumes  that  events 
on  a  link  are,  to  a  good  approximation,  independent  of  upstream 
offsets,  TRANSYT  carries  more  memory  of  previous  traffic  history. 

However,  TRANSYT' s  optimization  procedure  is  not  global  and  is 
computationally  rather  slow,  thereby  /limiting  the  number  of 
variables  that  can  reasonably  be  searched  over.  Finally,  we  note 
that  Gazis[12]  proposed  to  use  a  piecewise  linear  link  performance 
function  so  that  linear  programming  can  be  applied  to  calculate 
optimal  offsets  and  that  Eisner[6]  adopted  a  simple  "two-piece 
linear"  function  proposed  by  Newell[33]  and  formulated  the  network 
coordination  problem  in  terms  of  mixed-integer  linear  programming. 
However,  no  computations  have  been  attempted  by  these  two  workers. 
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On  the  hardware  side,  the  new  technology  has  produced  flexible, 
reliable  signal  switching  equipment.  A  complex  network  can  rather 
easily  by  controlled  from  a  central  location  with  a  small  computer. 
Such  a  computer  can  also  store  a  variety  of  pre-calculated  signal 
settinqs  and  put  them  into  effect  by  time  of  day  or.  to  a  limited 
degree,  in  response  to  inputs  from  vehicle  detectors.  The  new  systems 

have  the  important  characteristic  of  being  easy  to  change  [36]. 

A  new  timing  plan  can  be  introduced  by  reading  in  a  few  cards  at 
the  central  location. 

Real-time,  closed-loop  systems  that  employ  extensive  vehicle 
detection  and  powerful  central  computers  have  caught  the  imagination 
of  many  traffic  engineers  and  city  managers.  Such  systems  offer  the 
hope  of  making  the  best  use  of  a  severely  limited  resource,  namely, 
the  existing  city  streets.  A  variety  of  projects  have  been  under¬ 
taken,  starting  with  the  pioneer  work  in  Toronto  [14].  A  number  of 
worthwhile  benefits  have  been  obtained.  However,  such  systems  cur¬ 
rently  have  several  drawbacks,  including  high  Initial  cost,  con¬ 
siderable  ongoing  reliability  and  maintenance  problems  (particularly 
with  the  detectors)  and,  most  important  to  our  discussion,  a  lack  of 
adequate  theory  for  deciding  how  to  go  in  some  optimal  manner  from 
detector  input  to  signal  settings. 

The  present  work  is  part  of  a  research  and  development  program 
undertaken  by  the  Federal  Highway  Administration  in  connection  with 
the  UTCS  (urban  Traffic  Control  System)  in  Washington,  D.C.  The 
program  seeks  to  proceed  from  simple  systems  to  complex  ones,  through 
a  series  of  "generations."  Briefly,  first  generation  consists  of 
fixed-time  settings  calculated  offline  and  stored  in  simple  signal 
coordination  hardware.  Second  generation  also  uses  fixed-time  set¬ 
tings  but  calculates  them  online  in  response  to  traffic  detection. 
Third  generation  embodies  continuous  computer-calculated  control. 

The  research  approach  described  in  the  following  pages  has  been 
first  to  develop  a  traffic  model  and  performance  measure  and  then  to 
bring  contemporary  optimization  techniques  to  bear  on  the  problem. 

Tbe  principal  technique  used  is  mixed-integer  linear  programming 
and  so  the  method  has  been  given  the  acronym  MITROP  (mixed  integer 
traffic  optimization  program).  In  the  researchers1  vTew  the  results 
obtained  so  far  are  both  significant  and  promising.  The  emphasis  on 
a  systematic  approach  has  brought  new  advances  applicable  to  first 
and,  potentially,  second  and  third  generation  systems. 
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The  main  features  of  MITROP  are: 

•  the  simultaneous  consideration  of  all  major  control 
variables,  i.e.,  offsets,  splits,  and  cycle  time 

•  a  global  optimization  technique 

•  the  use  of  link  performance  functions,  and 

•  a  minimum  delay  objective  function  that  includes 
not  only  the  usual  deterministic  queue  buildup 
and  discharge  but  also  the  stochastic  effects  of 
overflow  queues.  The  latter  occur  as  flow  rates 
approach  intersection  capacity. 

The  characteristics  of  MITROP  can  be  brought  out  by  comparing 
it  to  the  Combination  Method  (CM)  and  TRANSYT.  Consider  Cfl  first. 
Both  CM  and  MITROP  use  link  performance  functions  but  MITROP 
permits  simultaneous  optimization  over  splits  and  cycle  time.  CM 
could  be  extended  to  handle  cycle  time  also  at  a  very  considerable 
increase  in  computation  time  (see  Gartner  and. Little  [11] 
but  the  inclusion  of  splits  is  virtually  out  of  the  question  (see 
Korsak[23]  ).  TRANSYT  contains  a  more  detailed  representation  of 
traffic  than  any  of  the  systems  mentioned  above.  However,  it  is 
net  yet  clear  how  significant  this  is  at  moderate  flows  since 
TRANSYT  and  SIGOP  have  been  found  comparable  in  field  evaluation 
both  in  Glasgow  and  San  Jose  [41  ,22]  .  Furthermore,  TRANSYT  uses 
one  variable  at  a  time  search  for  estimating  optimal  values  for  the 
decision  variables.  Convergence  is  to  a  local, not  global  optimum 
and  therefore  the  results  depend  on  the  starting  settings. 

Computation  is  relatively  slow  and  although,  in  principle,  all 
control  variables  (splits,  offsets,  and  cycle  time)  can  be  examined, 
in  practice,  search  has  been  confined  to  offsets,  primarily. 

MITROP  is  the  first  method  to  search  over  all  control  variables: 
offsets,  splits,  and  cylce  time  with  a  globally  optimal  method*  The 
most  useful  part  of  this  increased  scope  aopears  to  be  the  con¬ 
sideration  of  cycle  time.  Sensitivity  testing  on  sample  networks 
indicate  significant  improvements  in  performance  resulting  from 
.letting  cycle  time  vary.  For  example  the  optimal  settings  for  a 
particular  test  network  analyzed  in  the  report  when  all  variables  are 
considered  simultaneously  shows  an  improvement  of  11. 6X  in  overall 
performance  relative  to  settings  obtained  when  the  control  variables 
are  considered  separately. 
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.It  is  important  to  understand  the  physical  tradeoffs  involved  in* 
determining  cycle  times.  If  it  were  not  for  switching  losses,  very  short 
cycle  times  would  be  preferable.  This  is  because,  if  capacity  is 
adequate,  delay  is  never  more  than  one  red  period  and  decreases  with 
cycle  length.  However,  the  need  for  increased  capacity  favors  long 
cycle  times.  In  order  to  make  the  appropriate  tradeoff  on  a  network¬ 
wide  basis,  it  is  necessary  to  know  at  every  intersection  how  much 
worse  the  traffic  situation  is  when  cycle  time  is  shortened.  What 
happens  on  the  street  is  that  when  flow  is  too  close  to  capacity  , 
a  few  cars  will,  from  time  to  time,  fail  to  clear  the  intersection 
during  a  single  cycle  and  thereby  create  an  overflow  queue  a'long 
with  substantial  extra  delay.  Tf  flow  and  capacity  become  even 
closer,  overflow  queues  happen  more  frequently  and  delay  becomes 
worse.  MITROP  contains  a  representation  of  this  phenonmenon,  via 
the  saturation  deterrence  function,  and  thus  is  able  to  make  the 
tradeoffs  essential  to  calculating  optimal  cycle  times.  Recommended 
directions  for  future  research,  exploiting  the  MITROP  results,  are 
suggested  in  Section  8. 
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2.  THE  TRAFFIC  SIGNAL  NETWORK  TIMING  PROBLEM 


This  chapter  presents  the  traffic  signal  network  timing  problem  in  a 
general  form.  The  traffic  network  is  modeled  in  terms  of  nodes,  representing 
the  signalized  intersections,  and  links,  representing  sections  of  streets 
carrying  traffic  in  one  direction  between  two  intersections  (see  Fig.  2.1). 
The  variables  and  parameters  of  the  system  are  defined  and  the  physical  con¬ 
straining  relations  existing  among  them  are  stated.  The  timing  problem  is 
then  formulated  as  a  nonlinear  optimization  program  involving  the  minimiza¬ 
tion  of  a  generalized  disutility  function  associated  with  traveling  through 
the  signalized  intersections. 

2.1  Nomenclature 


The  traffic  network  and  the  principal  control  variables  are  defined 
as  follows: 


$„.  =  traffic  signal  at  node  i 

N  a  {i}  =  the  set  of  all  nodes  in  the  network 

(i,j)  =  the  link  going  from  node  i  to  j 

L  “  *  the  set  of  all  links  in  the  network 

^  =  [N,L]  =  the  network  graph 

Rij  =  Physical  rec*  time  at  facing  link  (i,j) 
r^j  =  effective  red  time  at  Sj  facing  (i,j) 

G-h  =  physical  green  time  at  S.  facing  (i,j) 

9ij  s  ®tfective  green  time  at  S..  facing  (i,j) 

Ajj  =  amber  time  at  facing  (i,j) 
ljj  *  phase  lost  time  at  Sj  in  direction  (i,j) 
a^  =  periods  of  all -red  time  at  Sj 
C  =  the  network  cycle  time 

All  time  variables  defined  above  are  measured  in  seconds. 

The  relations  between  physical  and  effective  signal  timings  are 
illustrated  in  Fig.  I. 2,  following  the  basic  model  for  traffic  signal  opera¬ 
tion  used  by  Webster  [40]  and  others  [5,27,39]. 
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Possible  flow 
through  the 
intersection 


GREEN 

AMBER 


G  +  A  *  g  ♦  l 
r  +  g  =  C 

1  *  V  +  1" 


r-j  -  (1,  +  a)  s  92  +  12 
r2  “  ^2  ■*■*)*  91  +  *1 


Figure  2.2  -  Basic  model  for  traffic  signal  operation. 
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Referring  to  Fig.  2.3  we  shall  now  relate  the  concepts  of  offset 
and  arrival  time  on  a  link.  Vehicle  platoons  released  during  the  green  time 
of  $•  travel  towards  S-.  Let 

■  J 


x- •  =  travel  time  of  the  platoon's  leading  edge  from  S.  to  S- 
J  on  the  connecting  link  1  J 

y.-i  =  arrival  time  of  the  platoon's  leading  edge  at  S.'s  step 

•  J  J 

line,  measured  relative  to  start  of  g..  so  that  for  every 
cycle  we  have,  J 


f*  •  ♦  ^  v .  •  ■<£  Q  *  • 


(2.1) 


offset  time  from  beginning  of  green  at  Sj  to  beginning  of 
respective  green  at  S.  so  that,  1 

J 


(2.2) 


=  offset  time  between  and  Sj  measured  from  beginning  of 

green  at  S.  to  beginning  of  next  green  at  S^.  Consequently, 
we  have,  1  J 


eij  ^ij-Wd  C 


0  _<  e.j  <  C 


(2.3) 


Finally,  let 


zij  ^ij*rij,9ij*C^  =  avera9e  cost  Per  vehicle  on  link  (i,j) 
when  traveling  through  S- 

fy  =  average  flow  on  link  (1,j)  (vehicles/second). 

2.2  Objective  Function 


is  a  disutility  incurred  by  vehicles  traveling  through  the 

signal -control led  network  such  as  delay,  stops,  acceleration  noise,  etc.,  and 
therefore  should  be  minimized.  Hence,  the  objective  function  in  an  optimiza¬ 
tion  procedure  would  be: 


min(i^),,Lfij  z1j  (*1j’rij’9ij»C) 


(2.4) 
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This  formulation  assumes  a  separability  condition,  i.e.,  the  dis¬ 
utility  on  each  link  is  a  function  of  parameters  for  that  link  only. 

Newell  [31]  has  shown  in  a  theoretical  study  that  this  condition  holds  tor 
high  traffic  flows  and  when  platoons  formed  by  the  traffic  lights  spread  so 
much  between  lights  that  the  time  between  the  arrival  of  the  first  and  last 
cars  of  a  platoon  at  a  light  exceeds  the  greer.  period.  The  validity  of  this 
assumption  is  also  supported  by  several  empirical  studies  [10, 16].  Further¬ 
more,  the  British  Combination  Method,  which  premises  the  same  assumption,  has 
shown  good  results  when  implemented  in  field  studies  [18]. 

The  constraints  which  the  optimization  procedure  must  obey  are 
described  next. 

2.3  Constraints 


Let 


i  =  a  set  of  links  that  form  a  loop  (undirected  circuit)  of  the 
network 

=  the  set  of  all  loops  of  the  network 
We  also  use  the  following  set-theoretic  notation: 
e  =  that  belongs  to 
V  =  for  every 

A  physical  constraint  of  the  system  is  that  the  sum  of  the  offsets 
around  any  loop  must  be  an  integer  multiple  of  the  cycle  time,  i.e., 

(A)**11 c  V1£<S£  (2-51 


where  n  is  an  integer  variable  associated  with  loop  t.  Another  physical 
constraint  is 


rij  +  9ij=C  V(i,j)eL  (2.6) 

In  order  for  the  network  to  be  able  to  handle  the  given  flow  we  have 
the  capacity  constraint: 

fij  Sij  -  fij  C  V  (i »j)  e  L  (2.7) 

where  s^j  is  the  saturation  flow  rate  at  the  downstream  signal  of  (i,j). 


To  facilitate  pedestrian  crossing,  minimum  red  must  be  established 

min 


ri j  -  rij 


(2.8) 
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Since  the  gain  in  capacity  with  vsry  long  cycle  times  is  often 
insignificant,  we  should  have  an  upper  limit  on  the  cycle  length.  This  is 
also  desirable  for  preventing  drivers  from  becoming  extremely  impatient  or 
believing  that  the  signals  are  defective. 


CiC1™"  (2,9) 

From  a  practical  point  of  view,  inching  safety  considerations,  it 
is  desirable  also  to  have  a  lower  limit  on  the  cycle  time  (i.e.,  an  upper 
limit  on  the  aspect  switching  rate),  in  addition  to  capacity  considerations, 

C  >  Cmin  (2.10) 


For  simplicity  at  this  stage  of  model  development,  we  mums  there 
are  no  special  turn  signals,  lagged  greens,  etc.,  and  the  intersections  in¬ 
volve,  at  most,  two  intersecting  streets.  However,  link  omissions  such  as 
one-way  streets  and  T  intersections  are  accommodated.  We  need  a  notation  to 
say  that  green  times  are  the  same  when  approaching  the  signal  from  eithc* 


direction  on  the  same  street  and  to  express  the  relation  between  red  on  one 
street  (e.g.,  the  main  street)  and  green  on  the  intersecting  street  (e.g., 
the  cross  street).  For  this  purpose,  let 


D 

j  «  {i  |  (i,j)  is  on  main  street) 


P-i  *  (1  |  (1  *d)  is  on  cross  street) 

J 


The  physical  red  and  green  times  are  -  14_-  and  g4.  +  14^,  respectively 
(Fig.  2.2).  10  1J 


The  equality  of  physical  reds  and  greens  on  opposing  approaches  to 
node  j  on  the  same  street  is  expressed  by 


r1j  “  11j  =  rk.i  "  1kj 

|  it* 

C-11) 

9ij  +  ]1j  =  9kj  +  1kj 

)  ■ 

For  the  equality  of  physical  red  and  green  on  a  street  and  Its  cross 

street  we  refer  to  Fig.  2.1  and  Fig.  2,2: 

r1j  -  <’lj  +  aj>  *  9mj  +  Vi  | 

I  j  e  N 

l 

V  *  <V  +  aj)  -  9lj  +  j 

it?, 
me  P4 

(2.12) 
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2.4  Optimization  Problem 

Summarizing  the  characterization  given  above,  the  network  signal 
setting  problem  can  be  stated  as  a  mixed- integer  nonlinear  optimization 
(MINLO)  program: 


MINLO:  Find  values  of  4^,  g^j,  C  to 
"1n(i?j)eLf1J  ZlJ  (*1J’  rij’  Slj’  C> 


subject  to: 

I  ♦ii  =  n#  C 
13  * 

rij  +  9ij  =  C 


VJt  e 


£ 


91j  sij  if1j  C>  v  c  L 


r  >  rmin 

r1j  -  rij 
rij  '  'ij  *  pkj 


-  1 


kj 


]  V  j  e  N 
\  V  i,k  e  Pj 


9ij  +  ’ij  *  9kj  +  ]kj  J 
rij  '  <’ij  *  aj>  9mi  +  ’mj  )  J  *1  l 
rmj  "  ^mj  +  aj J  *  9ij  +  ’ij  j  y  m  £  Pj 

(jtiin  <  j.  <  (jnax 


♦ij,n£  unrestricted  in  sign;  n£  integer. 

A  substantial  number  of  the  variables  and  constraints  are  redundant  and  can 
be  eliminated  by  substitution. 

Among  the  various  candidate  optimization  techniques  for  solving  this 
program  we  have  chosen  to  use  mixed-integer  programming  (KIP).  One  of  the 
chief  attractions  of  this  approach  is  that  reasonably  efficient  packages  for 
solving  such  problems  have  recently  become  available,  especially  MPSX  which 
runs  on  IBM  370  equipment  [21]. 
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3.  LINK  PERFORMANCE  FUNCTION 


The  disutility  associated  with  traveling  through  the  signalized  links  is 
described  in  this  section  in  terms  of  a  Link  Performance  Function  (LPF).  A 
traffic  flow  model  is  developed  for  this  purpose  and  evaluated  with  respect 
to  field  data. 

3.1  Basic  Relations 


In  the  following  discussion  we  establish  the  beginning  of  green  time 
at  the  downstream  signal  of  the  link  as  a  reference  point  and  assign  to  it 
zero  value.  Thus,  the  time  interval  [-r,g]  denotes  one  cycle  period  con¬ 
sisting  of  an  effective  red  period  [-r,0]  and  an  effective  green  period 
[0,g]  (Fig.  3.1).  This  section  deals  exclusively  with  traffic  behavior  on 
a  specific  link  (i,j).  To  simplify  notation  we  drop  subscripts. 


i 

one  cycle 

i _ _ _  j _ 

i _ l _ 

1 

i  i  :  i  •  n 

0  .9  time 

Figure  3.1:  Signal  parameters. 


Definitions: 

q(t)  =  arrival  rate  (vehicles/second) 

A(t)  =  cumulative  number  of  arrivals  at  time  t  (vehicles). 

Starting  at  the  beginning  of  any  red  period,  we  have 

A(t)  =  / 1  q(r)  dv  (3,1) 

-r 

s  =  saturation  flow  rate  at  the  signal's  stop  line  (vehicles/ 
second) 

s(t)  =  possible  departure  rate  at  time  t.  For  each  cycle  we  have: 

!0  for  -  r  <  t  <  0 

(3.2) 

s  for  0  <.  t  <  g 

The  saturation  flow  s  is  maintained  during  the  green  period  only  as  long 
as  there  is  a  queue  waiting  for  service.  Otherwise,  the  departure  rate  equals 
the  arrival  rate. 
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We  first  confine  ourselves  to  deterministic  flows.  Stochastic  effects 
are  accounted  for  later  in  Section  5.  We  assume: 

1.  Arrivals  are  periodic,  i.e., 

q(t)  =  q(t  +  n  C)  (3.3) 

where  n  is  an  integer  number. 

2.  The  signal  is  undersaturated,  i.e., 

A_  <  gs  (3.4) 

where  v 

Ad  =  I9  dt  (3.5) 

F  *r 

is  the  total  number  of  cars  arriving  at  the  signal  during 
one  cycle,  ys  denotes  the  capacity  of  the  signal,  that  is, 
the  number  of  possible  departures  during  the  green  period 
(and  hence  also  during  one  cycle  time). 

3.  The  arrival  rate  during  the  green  period  of  the  signal 's 
cycle  time  does  not  exceed  the  saturation  flow  rate,  i.e., 

q(t)  <  s  for  0  <  t  <  g.  (3.6) 

(3.6)  implies  that  once  a  queue  has  vanished  during  the 
green  period  it  cannot  rebuild  before  the  next  red  period 
commences. 

According  to  these  assumptions,  all  vehicles  arriving  during  a  cycle 
in  which  the  red  period  precedes  the  green  can  be  accommodated  in  that  cycle. 

It  follows  that  the  queue  is  always  empty  at  the  end  of  the  green  period  and 
delay  time  calculations  can  be  confined  to  a  single  interval  [-r,gj. 

The  queue  length  Q(t)  at  any  time  t  e  [-r,g]  is  given  by  the  dif¬ 
ference  between  the  cumulative  number  of  arrivals  and  the  cumulative  number 


departures: 

(  A(t) 

if  -r  £  t  <  0 

Q(t)  =  <  A(t)  -  ts 

if  0  £  t  £  tQ 

(3.7) 

(° 

if  t  <  t  <  g. 
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tQ  is  the  queue  clearance  time  and  denotes  the  first  time  during  the  green 
period  when  the  queue  disappears.  By  definition  t  =  tQ  when 

Q(t)  =  A(tQ)  -  tQ  s  =  0  and  0  <.  tQ  <.  g.  (3.8) 

The  delay  incurred  by  Q(t)  queueing  vehicles  during  an  interval  dt  is 
Q(t)  dt.  Hence,  the  total  delay  time  2  incurred  by  traffic  during  one  cycle 
is  represented  by  the  area  under  the  queue  length  curve: 

Z(Y,r)  =  J9  Q(t)  dt  =  Z0  Q(t)  dt  (3.9) 

-r  -r 

The  average  delay  per  vehicle  z(ytr)  is  obtained  by  dividing  by  the  total 
number  of  arrivals  during  one  cycle: 

z(Y,r)  =  J-  Z(y*r)  (3.10) 

AP 

Note  that  z  is  a  function  of  arrival  time  and  split  and  is  independent  of 
cycle  time. 

The  subsequent  discussion  considers  primarily  delays  as  the  cost  to 
be  minimized.  However,  the  models  that  are  developed  can  be  used  also  to 
calculate  a  more  general  link  performance  function  combining  costs  of  delays, 
stops,  acceleration  noise,  or  other  measures  of  effectiveness,  by  using 
appropriate  weighting  factors  [4,19]. 

3.2  Traffic  Flow  Model 


We  construct  a  simple  model  to  analyze  the  link  performance  function. 
We  assume  the  traffic  flow  pattern  on  the  link  to  be  a  rectangular-shaped 
periodic  platoon.  The  platoon  is  characterized  by  its  time-length  p  and 
constant  flow  rate  q  (Fig.  3.2). 


P 


Figure  3.2:  Traffic  flow  pattern. 
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We  can  always  determine  for  this  model  a  period  of  one  cycle  time  In  which: 


(q  if  Y<t<y  +  p 

q(t)  =  (3.11) 

(0  else. 

The  platoon  size  for  this  model  is  given  according  to  (3.5)  by: 

Ap  *  qp  (3.12) 

Referring  to  Fig.  2.1  we  describe  how  to  determine  the  traffic  flow 

parameters  on  link  (j,k).  Given  the  effective  green  time  at  the  upstream  end 
of  link  (j,k),  g,.,.»  we  assume  that  the  platoon  length  p<u  when  entering  link 
(j,k)  Is:  1J  Jlt 

Pjk  *  g,.j.  (3.13) 

If  the  average  flow  on  the  link  is  f4k  (vehicles/second),  we  further  assume 
that 


qjk 


fjk  9lj 


(veh/sec). 


(3.14) 


i.e.,  the  flow  is  assumed  to  be  uniformly  distributed  over  the  entire  green 
phase.  Alternative  assumptions,  such  as  a  'tadpoling*  flow  pattern  [20],  or 
other  patterns,  could  be  equally  made. 

Equations  (3.13)  and  (3.14)  refer  to  the  primary  flow  on  link  (j,k). 
If  there  is  a  substantial  secondary  flow,  the  platoon  length  is  modified  as 
shown  in  Fig,  3.3  (we  use  now  superscripts  1  and  "  for  primary  and  secondary 
flows,  respectively).  Consequently  we  have, 

p  =  p'  +  p".  (3.15) 

p'  is  determined  solely  by  the  green  phase  serving  the  primary  flow  on  link 
(j,k),  i.e.,  pjk  =  g^.  as  in  eqn.  (3.13)  above.  Similary,  qjk  Is  determined 

following  (3.14), 


=  fjk  <veh/sec> 


(3.16) 


Figure  3.3: 


Primary  and  secondary  flows. 


In  order  to  maintain  a  uniform  platoon,  the  extended  portion  due  to  the  second¬ 
ary  flow  is  calculated  as  follows: 


Hence, 


(3.17) 

(3.18) 


In  certain  cases  the  turning-in  traffic  constitutes  the  primary  flow,  for 
instance,  the  flow  from  (m,j)  onto  (j,k)  in  Fig.  2.1  may  be  the  primary  flow 
on  link  (j,k).  In  these  cases  p'  is  determined  by  the  appropriate  green  time 
from  which  it  emanated  (g  -  in  this  example). 


To  determine  the  link  performance  function  we  need  p  and  q  at  the 
downstream  end  of  the  link.  Both  parameters  are  a  function  of  the  link 
length  d,  owing  to  the  dispersion  phenomenon.  A  simple  relationship,  cor¬ 
roborated  by  some  empirical  evidence  [13,30]  is  assumed: 

p(0)  x  k(d)  if  p(d)  <  C 

P(d)  =  (3.19) 

C  if  p(d)  >  C 
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where  k(d)  is  a  piecewise  linear  relationship  as  shown  in  Fig.  3.4.  The 
platoon  length  cannot  exceed  the  length  of  one  cycle  time,  at  which  point 
flow  is  becoming  continuous  (see  also  [42]). 


k(d) 


d ( ft . ) 


Figure  3.4:  Platoon  dispersion  factor. 


Relationships  that  are  different  from  the  one  shown  in  Fig.  3.4  may  be  taken 
on  specific  links,  if  manifested  so  by  empirical  evidence  (one  different 
example  is  demonstrated  by  the  experimental  data  shown  in  the  next  section). 
q(d)  is  determined  by  flow  conservation  considerations: 


therefore. 


q(d)  p(d)  =  q(0)  p(0), 

q(d)  -  q<0)  {$j}  • 


3.3  Evaluation  Results 


(3.20) 


The  traffic  flow  model  described  in  Section  3.2  was  evaluated  by  com¬ 
parison  with  field  data  on  platooned  traffic  reported  by  the  British  Transport 
and  Road  Research  Laboratory  [16].  In  this  experiment,  observations  on  in¬ 
dividual  vehicle  arrivals  were  recorded  at  four  different  positions  downstream 
of  a  signal -controlled  intersection  operating  on  a  90-second  cycle.  This  in¬ 
formation  was  subsequently  processed  to  determine  the  average  number  of  arrivals 
in  each  two-second  time  interval ,  referred  with  respect  to  the  onset  of  the 
green  phase  at  the  upstream  traffic  signal,  for  each  one  of  the  four  locations. 
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The  results  together  with  a  scaled  schematic  of  the  intersection  and  roadway 
where  the  information  was  collected  are  shown  in  Fig.  3.5.  The  spreading, 
or  dispersion,  that  takes  place  throughout  the  approximate  1000  ft.  of  roadway 
downstream  of  the  intersection  can  be  readily  seen,-  e.g.,  the  fourth  platoon 
has  72-sec.  passage  time  that  should  be  compared  to  the  40-sec.  effective 
green  time  at  the  intersection. 

The  actual  platoons  were  approximated  by  rectangular-shaped  platoons 
containing  the  same  total  number  of  vehicles.  The  width  was  made  correspondent 
to  a  pair  of  imaginary  rays  delineating  the  trajectories  of  the  leading  edge 
and  the  trailing  edge  of  the  rectangular  approximation.  The  relationships  of 
Section  3.1  were  then  used  to  calculate  average  delay  per  vehicle  in  seconds 
z(y)  as  a  function  of  arrival  time  y,  parameterized  for  three  splits  g/r  for 
both  the  actual  platoons  and  their  rectangular  approximations.  The  splits 
were  chosen  to  cover  a  wide  range  of  possible  degrees  of  saturation.  The 
results  are  illustrated  in  Figs.  3.6  to  3.9,  corresponding  to  each  of  the 
four  different  platoons  of  Fig.  3.5.  The  following  observations  can  be  made: 

In  most  cases  there  exists  a  close  fit  between  the  delay  functions  derived 
from  the  actual  platoons  and  those  derived  from  their  rectangular  approxima¬ 
tions.  The  largest  deviations  occur  at  extremes  of  splits  and/or  arrival  times 
(and  hence  offsets)  that  are  only  rarely  expected  in  practice.  In  some  of 
these  cases  it  is  possible  to  improve  the  fit  through  a  better  selection  of 
the  parameters  y  and  p  for  the  rectangular  platoon.  Clearly,  further  exper¬ 
ience  and  field  data  ought  to  guide  this  selection. 

3.4  Analysis  of  the  LPF  for  Rectangular  Platoons 

In  this  section  further  relations  are  developed  to  characterize  the 
link  performance  function  z(y,r),  based  on  the  particular  traffic  flow  model 
described  above.  These  relations  will  be  useful  in  developing  the  analytics 
to  be  employed  later  in  the  network  optimization  program.  As  was  mentioned 
previously,  the  approach  outlined  in  this  section  can  be  similarly  applied  to 
additional  measures-of-effectiveness,  such  as  number  of  stops,  acceleration 
noise,  etc. 

For  a  fixed  platoon  length  p,  the  average  delay  z  is  a  function  of 
arrival  time  y  (translatable  to  offset  time  <j>  via  equation  (2.2) )  and  the 
split  (say  r).  We  can  prove  the  following  theorems: 

Theorem  3.1:  Delay  is  minimal  when  the  arrival  time  y  is  such  that  the 
trailing  edge  of  the  platoon  arrives  at  the  light  just  before  it  turns  red 
and  is  able  to  clear  the  intersection  entirely,  i.e., 

y  +  p  =  g.  (3.21) 

Proof:  Referring  to  Fig.  3.10  we  distinguish  two  cases: 

(a)  p^g,  (b)  p  >  g. 
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Figure  3.5:  Average  platoon  flow  and  rectangular 
approximation  at  four  observation 
posts  downstream  of  a  signal -controlled 
intersection. 
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Case  a:  A  degeneracy  occurs  for  y  e  [0,g-p].  The  platoon  fits  into  the 
green  opening  and  since  no  cars  are  stopped  the  delay  is  zero 
(case  a.l).  For  y  <  0  we  have  a  headstoppinq  effect:  q(-y) 
cars  arrive  during  red  and  the  delay  incurred  is  indicated  by 
the  dashed  area  shown  in  case  a. 2.  For  y  +  p  >  g,  a  portion 
p„  of  the  platoon  is  tail  chopped  by  the  red  light:  qp„  cars 

A  A 

are  stopped  and  delayed  (case  a. 3). 

Case  b:  When  y  +  p  =  g,  i.e.,  y  =  -(p-g),  q(-y)  are  stopped.  The  delay 

is  indicated  by  the  dashed  area  shown  in  case  b.l.  For  y  <  -(p-g) 
headstopping  increases  and  delay  increases  correspondingly 
(case  b.2).  For  y  >  -(p-g)  the  total  number  of  cars  that  arrive 
during  red  is  the  same  as  in  case  b.l,  i.e.,  q(p-g).  However, 
in  this  case  both  headstopping  and  tail  chopping  occur:  qpx  are 

stopped  right  after  the  light  turns  red  and  thus  the  total  delay 
is  larger  than  in  b.l  (case  b.3).  Q.E.D. 

Referring  to  Fig.  3.10  (case  b.2)  we  calculate  the  total  delay  Z(y) 
incurred  by  the  platoon.  It  is  the  sum  of  two  areas: 

£ 

Z(y)  =  A-j  +  A2  =  \  qy2  +  /  °[q(t-y)  -  st]dt. 


The  queue  clearance  time  tQ  is  determined  by  the  identity 

Q(t0)  =  q(-y)  +  qtQ  =  stQ  =  0, 
following  eqn.  (3,8).  Therefore, 


t  =  _91  =  _Y  JL. 
o  q-s  Y  1-y’ 


(3.22) 


where  y  =  q/s  denotes  the  ratio  of  platoon  flow  to  saturation  flow  on  a 
particular  link.  Hence, 

a2  =  T(q  -  s)  I  -  qyt]0°  =  - 


and. 


Z(Y)=4or2-!-s^?  =  ^a§^ 


2  qr 

The  average  delay  is  given  by 

Z/Y)  =  Ihl  =  A  .  _L 

qp  2p  1-y' 


(3.23) 
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By  Theorem 
Therefore, 


3.1  the  minimal  value  occurs  when  y  =  x  = 


zmin 


1 

1-y 


1 

i-y 


-(p-g). 


(3.24) 


Theorem  3.2:  Delay  is  maximal  when  the  arrival  time  y  is  such  that  the 
leading  edge  of  the  platoon  arrives  at  the  signal  stop  line  right  after  the 
light  turns  red,  i.e.. 


y  =  y  =  -r  (3.25) 

Proof:  Referring  to  Fig.  -3.11  we  distinguish  three  cases: 

(a)  p  <  r,  (b)  r  <  p  <  tQ  +  r,  (c)  p  >  tQ  +  r. 

Case  a:  Fig.  a.l  shows  the  queue  evolvement  and  total  delay  for  y  =  -r. 

For  y  >  -r  (case  a. 2)  the  platoon  arrives  later  during  the  red 
period  and  the  stopped  cars  wait  a  shorter  time  for  the  green 
light  than  in  case  a.l.  Similarly,  for  y  >  0,  a  portion  (g-y) 
of  the  platoon  clears  the  intersection  before  red  starts  and 
delay  is  reduced  in  this  case  as  well  (case  a. 3). 

The  total  delay  in  case  a  can  be  described  as  the  summation  of  three 
sub-areas: 

A1  =  ?  qp2;  A2  =  -W(r*l>)s  a3  =  /  °(pq  -  st)dt. 

The  clearance  time  tc  is  determined  by  the  identity  pq  =  stQ;  i.e., 

t0  =  fp  =  yp  (3.26) 

1  2 

Consequently,  A3  =  ^qyp  and  the  total  delay  is, 

Z(y)  =  A1  +  A2  +  A3  =  qp(f--  |  -  y). 

Following  eqn.  (3.10), 

z(a)  (y)  =  f(y-l)  -  y  (3.27) 

indicating  that  the  average  delay  decreases  linearly  as  y  is  increased.  The 
maximum  delay  occurs  at  y  =  -r. 


Figure  3.11-(a) 


I 

Pr 
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(3.28) 


Case  b:  Analysis  is  similar  to  that  of  case  a(see  Fig.  3.11  -  (b) ) . 

The  maximal  total  delay  in  this  case  is  again  a  summation  of  three 
sub-areas*  where 

A,  =  /  q(t-y)dt  = 

Y 


P+Y  ^2 

A?  =  /  [q  (t-Y)  -  St]dt  =  (q-s)  -  qy(p+Y) 


A3  =  /  °[Qt=p+Y  -  s(t-P-Y)]dt. 

pH 

The  queue  length  at  t  =  p+y  is  Q+_n.  =  qp  -  s(p+y)  and  tn  =  yp  as  in  case  a. 
Consequently,  0 

2 

A3  *  \  -  qp(p+Y)  +  ]s(p+y)2. 


By  addition  we  obtain. 


7,  x  _  1  (qpr  ,2  12 

Z(y)  =  j  A  y  ■  ■  -  ^(p-hr)  +  t<1y 


z(b)(T)  =  lp[yp2  -  (p+y)2  +  Y2]  *  f(y-D  -  T  -  2(al(T). 

Therefore,  for  y  =  t  =  -r  we  have 

,«•) .  -  z 

max  max  max 


(3.29) 


(3.30) 


(3.31) 


Case  c:  p  >  tQ  +  r,  i.e.,  flow  is  such  that  arrivals  extend  beyond  the 
queue  clearance  time.  Substituting  tQ  =  yp  we  obtain 

P  >  (3.32) 

(The  same  condition  is  obtained  when  substituting  expression  (3.33)  derived 
for  tQ  below.) 
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Analysis  is  again  similar  to  cases  (a)  and  (b).  Referring  to 
Fig.  3.11  -  (c),  total  delay  is  a  sum  of  two  areas: 


1-  2 


A-j  =  /  q(t-Y )dt  =  ^qy 


A2  =  /  °[q(t-y)  “  st]dt  =  ^qY2  • 


tQ  is  determined  according  to  eqn.  (3. ft): 

qy  +  qtQ  -  tQS  =  0 


t  =  -iLy  =  JLV 
''o  s-q’  1  -y Y  * 


(3.33) 


By  addition  we  obtain. 


and. 


Z(y)  *  +  a2  *  JjflY2 

2  1 

z(y)  =  in  ‘  Ty- 


1 


2p  i-y 

To  obtain  the  maximal  average  delay  for  this  case,  we  substitute  y  =  -r: 


(3.34) 

(3.35) 


*(«>  =  !_.  i 

nnax  2p  1-y 


(3.36) 

Q.E.D. 


Referring  to  Fig.  3.11-(c.3)  we  observe  that  a  degenerate  case  occurs 
for  pv  >  r/(l-y),  where  pv  *  y+p-g.  zill  maintains  a  constant  value  in  this 

A  A  liiaA 

case  for  as  long  as 


The  different  cases  illustrated  in  Fig.  3.10  and  Fig.  3.11  do  not 
cover  all  the  possible  combinations  of  arrival  time,  platoon  length,  split 
value,  flow  and  saturation  flow.  Two  generalized  analytical  expressions  for 
delays  incurred  by  platoons  at  traffic  signals  are  given  in  Ref.  [2]. 

We  are  now  able  to  deduce  the  following  two  corollaries,  confining 
ourselves  to  the  span  of  one  cycle  time,  i.e.,  y  e  [-r,+g]  (see  Fig.  3.12): 
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Corollary  1:  The  locus  of  minimal  delay  in  plane  (ytr)  is  the  line 
y  =  g-p  =  C-r-p.  Also,  delay  is  minimal  (and  equal  to  zero)  in  the  area 
circumscribed  by  0  £  y  <_C-r-p. 

Corollary  2:  The  locus  of  maximal  delay  in  plane  (y,r)  are  the  two  lines 
y  =  -r  and  y  =  g  =  C-r,  since, 

W  *  *<*•*•>  -  W.D.  '  (3-38) 

As  shown  above,  a  degeneracy  occurs  when  p  >  r(l-y).  For  any 
y  e  [g  -  (p-(r+tQ)),g]  there  is  no  change  in  delay  (at  a  fixed  split) and  its 

value  is  maximal  according  to  equation  (3.31). 

Theorem  3.3:  For  y  e  [-r,0]  the  delay  function  z(y,r)  is  a  .strictly  decreasing 
function  of  y. 

Proof:  Referring  to  Fig.  3.11  (cases  a.2,b.2  and  c.2),  any  positive  shift 
in  y  from  y^  to  y2»  such  that  y2  >  y-j»  results  in  a  queue-length  curve  Qg(t) 

that  is  completely  encompassed  by  Q-j(t).  Therefore,  z(y2,r)  <  z(y-j,r). 


Theorem  3.4:  For  y  e  [0,g]  the  delay  function  z(y,r)  is  a  non-decreasing 
function  of  y. 


Proof:  We  distinguish  the  following  cases: 


<c/(.1-y) 


>g 

<g 


i 

ii 


hi 

IV 


I.  Referring  to  Fig.  3.10  (cases  a. 3  and  b.3),  any  positive  shift  in 
y  from  y-j  to  Y2»  such  that  y2  >  y-j,  results  in  a  queue-length  curve 

Q2(t)  that  completely  encompasses  Q^(t).  Therefore,  z(y2,r)  >  z(y^,r) 

and  z(y,r)  is  strictly  increasing. 


II.  z(y,r)  has  a  zero-valued  flat  portion  for  y  e  [0,g-p]  and  is  strictly 
increasing  for  y  e  [g-p,g]. 


III.  z(y,r)  is  strictly  increasing  for  y  e  [0,C-p*ry/(l-y)]  and  has  a 

maximal  flat  portion  (variable  with  split  but  fixed  at  any  fixed  split) 
for  y  e  [C-p+ry/(l-y),g]. 

IV.  z(y,r)  is  strictly  increasing  for  y  c  Cg-p,C-p+ry/(l-y)]  and  has  two 
flat  portions  as  in  (II)  and  (III)  above. 
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Theorem  3.5:  If  p  >  r/(l-y)  then  for  y  g  [y.yj  delay  is  independent  of  split. 

Proof:  The  theorem  is  proved  most  easily  with  reference  to  Fig.  3.13.  The 
total  delay  incurred  by  the  platoon  in  cases  (a)  and  (b)  depends  only  on  the 
characteristics  of  the  platoon  itself  (p  and  g)  and  its  arrival  time  y. 

Since  all  these  parameters  are  identical  in  both  cases  (a)  and  (b),  the 
resultant  queue-length  curve,  and  hence  the  total  delay,  are  also  identical, 
regardless  of  the  split  value  (provided,  of  course,  we  do  not  exceed  capacity). 

Special  Case:  y  >  1 

The  discussion  in  this  section  assumed  up  to  now  that  y  <  1,  i.e., 
q  <  s.  Though  rare,  a  situation  may  occur  in  which  y  >  1  on  a  particular 
link.  For  instance,  in  case  a  two-lane  link  is  feeding  traffic  onto  a  single¬ 
lane  link  (the  bottleneck  occurs  at  the  intersection  connecting  them). 

The  capacity  constraint  requires  as  before: 

qp  <_  gs  i.e.  g>yp. 

Clearly,  since  y  >  1,  g  must  be  larger  than  d. 

(Caution:  Condition  (3.32)  on  p.41  now  has  to  be  read  p(l-y)  >  r,  since 
dividing  by  1-y  which  is  now  negative  would  change  the  sense  of  inequality.) 

Maximal  delay  is  given  by  equations  (3.28)  or  (3.36)  as  above. 

However,  calculation  of  minimal  delay  is  different.  Referring  to  Fig.  3.14 
we  have  a  queue  build-up  even  though  arrivals  occur  during  green  only.  Total 
delay  is  a  sum  of  the  two  areas  Aj  and  A2: 


Y+p  2 

A1  =  /  (q-s)(t-y)dt  =  (q-s)| 
y 

t  2 

Ag  =  /  [(q-s)p  -  s(t-y-p)]dt  =  (q-s)f  (y-1) 
Y+P 


where  tQ  is  determined  as  follows: 

q  •  P  -  s(t0y)  =  0 
=  Y+yp 

Consequently, 


4m =  z"2x(q-s> 


(3.39) 


(3.40) 
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4.  THE  TRAFFIC  SIGNAL  NETWORK  COORDINATION  PROBLEM 


In  this  chapter  we  formulate  the  traffic  signal  network  coordination 
problem  as  follows:  given  a  common  cycle  time  for  the  network  and  red  and 
green  splits  at  each  node,  find  a  fundamental  set  of  offsets  to  minimi2e 
total  delay.  The  formulation  is  in  terms  of  mixed-integer  linear  program- 
ming  and  is  followed  by  an  example. 

4.1  Piecewise  Linearization  of  the  Link  Performance  Function 

The  model  developed  in  Section  3  exhibits  delay  as  a  nonlinear 
function  of  offset,  split  and  cycle  time.  However,  if  delay  is  a  function 
of  offset  only  (split  and  cycle  time  are  kept  fixed  in  the  present  formula¬ 
tion)  we  have  a  quasi-convex  relationship  (see  Figs.  3. 6-3. 9).  To  accom¬ 
modate  mixed-integer  linear  programming  computational  techniques  it  is 
essential  to  linearize  the  objective  function  (total  delay  in  the  network) 
by  piecewise  linear  approximations.  Our  convention  is  to  use  secant 
approximating  lines  rather  than  tangent  lines.  They  seem  to  give  a  fairly 
accurate  approximation,  particularly  in  the  region  of  the  minimum,  and 
they  are  easy  to  construct.  Two  examples  of  such  an  approximation  of  delay 
versus  offsets,  for  two  different  splits,  are  shown  in  Fig.  4.1  and  Fig.  4.2. 

We  shall  now  transfer  the  linearized  decision  function  into  the 
constraint  set.  Suppose  we  approximate  the  objective  function  by  n  linear 
segments.  We  represent  each  line  as  z  =  z-,(<j>)  where  Z-.U)  is  linear  in  4. 

Now  consider  the  set  of  constraints 

z  >  Zj(*) 

z-z2^  (4.1) 

* 

z  1  zn(*)- 

They  imply  that  z  belongs  to  a  closed  convex  set  .  In  our  formulation 
we  introduce  one  linear  constraint  per  link  per  linear  segment  which  approxi¬ 
mates  the  true  objective  function.  This  accounts  for  the  bulk  of  the  con¬ 
straints  of  the  optimization  problem.  The  number  of  segments  depends  upon 
the  desired  accuracy  of  the  approximation.  Yet  one  should  keep  in  mind 
that  increasing  the  number  of  segments  increases  the  computational  burden. 

Another  property  of  delay  versus  offset  should  be  noted.  Suppose 
we  have  a  certain  delay  z-j  for  an  offset  $-j.  If  we  change  this  offset 

by  ?n  integral  multiple  of  the  cycle  time  we  should  obtain  the  same  delay, 
since  delay  as  a  function  of  offset  is  periodic  with  period  C.  Yet  we  have 
modeled  delay  as  a  nonperiodic  function  to  maintain  convexity.  Although 
the  function  is  not  explicitly  periodic,  the  integer  variables  in  the  offset 
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loop  constraints  reflect  this  property.  Suppose  we  have  the  piecewise 
linear  convex  objective  function  illustrated  in  Fig.  4.3.  Note  that 
z(x-g)  =  z{ t+r).  We  can  prove  the  following  theorem: 


i(<P) 


Figure  4.3:  Piecewise  Linear  Convex  Objective  Function. 


Theorem  41;  Any  optimal  value  of  offset  4>  must  belong  to  the  interval 
[x-g,  x+rj. 

Proof:  Suppose  not.  Let  z*  =  z(«j>*)  be  the  associated  delay,  with 

T+r].  There  exists  an  integer  m  such  that  <{>*  +  m  C  e  [T-g,  r+r] 
and  z(<j>*  +  m  C)  <  z since,  z(<l>*)  >.  z(x+r)  =  z(x-g)  >_  z(<^^  +  mC).  Further¬ 
more,  the  only  physical  constraint  in  which  offset  appears  is  the  loop  con¬ 
straint,  which  is  not  violated  by  changing  any  link  offset  by  an  integer 
multiple  of  cycle  time.  Thus  we  maintain  feasibility  and  consistency  of  all 
physical  constraints  and  the  optimization  procedure  must  produce  a  value  of 
z*  =  z($*)  such  that  <f>*  e  [x-g,  x+r].  Q.E.D. 

The  next  section  considers  in  more  detail  the  network  loop  con¬ 
straints. 

4.2  The  Network  Loop  Equations 

The  signal  network  loop  constraints  were  stated  in  a  general  form 

in  (2.5): 

l  *ij  =  ntC  v  1  e  (4.2) 

(i »j)  e  i 
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Let  us  consider  first  loops  containing  two  links  in  opposing  direc¬ 
tions  between  the  same  pair  of  nodes.  Loops  containing  two  links  in  the 
same  direction  are  a  trivial  case.  Referring  to  Fig.  2.3  we  obtain  the 
following  constraining  equation: 


<J>ij  $j.j  =  n  C  (4.2) 

Dividing  through  by  C  we  obtain: 

*ij  +  =  n  (4.3) 


where  <j>'  =  £  is  a  dimensionless  offset  variable  (given  in  fractions  of  cycle 

time).  The  integer  n  reflects  a  periodicity  in  the  variable  $  (or  $'), 
with  period  C.  Since  the  link  performance  function  as  a  function  of  offset 
has  also  period  C,  we  may  restrict  offsets  to  an  interval  of  length  C  and 
obtain  bounds  on  the  permissible  values  for  n. 


x . .  —  a . .  <  a . .  <  x-  -  +  rj* 

U  -  U  ij 

Thus,  taking  the  lower  and  upper  limits,  we  get  according  to  (4.2) 
<Tij  -  Sij>  +  <Tj1  -  9j1>  i  "  C  i  f'lj  +  r1j>  +  <Tj1  *  rj1> 


Defining, 
and  rearranging  (4.5)  we  obtain 


T1i  ”  Tij  +  Tj1 


Tjj  <  "  C  <  t1;j  +  2C 


or 

(t^ j/C)  1  n  ^  (t^/C)  +  2 

Let. 


Since 


and. 


[x]  =  {the  greatest  integer  less  than  or  equal  to  x} 
n  is  restricted  to  the  set  of  integers,  we  have, 
[t^/C]  +  1  i  n  <_  [t^/C]  +  2  if  L/C  is  not  integer, 

C*ij/C]  <  n  <  [i^/C]  +  2  if  j/C  is  integer. 

If  C  is  allowed  to  be  variable  as  in  Section  6  below 


(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 
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(4..12) 


cm1n  <  c  <  cma\ 

we  obtain  for  (4.10)  and  (4.11)  respectively: 

t^j/c”3*]  +  1  <  n  <  [;ij/Cnin]  +  2  (4.13) 

and, 

t^/C1"3*]  <  n  <  Pij/Cr“,n]  +  2  (4.14) 


For  the  case  where  a  loop  contains  more  than  two  links  we  mi:st  ex¬ 
tend  the  notion  of  offset.  So  far  ^ ■  has  been  defined  as  the  time  measured 

from  start  of  green  at  S-  to  start  of  green  at  Si  along  a  link  (i,j)  con- 

necting  the  two  signals.  That  is,  the  offset  is  measured  between  two  dif¬ 
ferent  nodes  in  the  same  direction  (see  Fig.  2.1).  Now  we  shall  also  con¬ 
sider  an  offset  representing  the  difference  in  the  start  of  green  times  at 
the  same  node  in  two  different  directions.  This  we  shall  term  an  intra-node 
offset;  the  time  representing  the  difference  between  the  start  of  green  at 
Sj  in  direction  (i,j)  and  the  start  of  green  at  in  direction  (j,m).  This 

difference  is  actually  the  effective  green  plus  lost  time  at  in  direction 
(i»j)»  or,  g^.  +  1-j  (Fig.  2.2).  Distinguishing  between  link  offsets  (inter¬ 
node  cffsets)  and  intra-node  offsets,  the  constraints  (2.5)  can  be  rewritten 
as 


,.^•3  (♦ii  +  'J'i-i)  =  n.  C  He  ^  (4.15) 

(i,j)  e  £  u  x 

where  =  g^.  +  1-  in  this  case. 

According  to  (2.6)  we  can  substitute  the  relation  g...  =  C  -  r._.  to 

obtain 


ah)  1 1  +  V  ’  (n* ' m)  c 

where  m  is  the  number  of  arcs  in  this  particular  loop  (m  >  2). 
Letting  n£  -  m  =  n  (integer)  we  obtain. 


(4.16) 


l 

(i»j)  e  £  V9ij  "  rij  T  ' i j 


•  +  l*,)  -  n  C 


(4.17) 
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and,  in  dimensionless  quantities, 

l  ($<,-  -  ri ,  +  w  lii)  =  n 
(i,j)  c  A  13  13  13 


(4.18) 


,  -1. 


where  w  =  1/C  denotes  the  signal's  frequency  [time  ]. 


To  estimate  the  integer  bounds  we  substitute  in  (4.16)  the  lower 
and  upper  limits  in  <)>,  as  in  (4.4)  above: 


<4  .  ^  +  ±nC~  (i^j)  e  +  '  r«  + 

¥  IeSO  ( 


Rearranging  and  substituting 


T,  -  .1 


*  (i.J)  e  £ 

we  obtain  for  the  bounds, 

v. 

y  l  e  o 
or , 

(yc)  -  m  <  n  <  (yc) 

For  C  variable  as  in  (4.12),  we  obtain, 

(Tt/Cmax)  -  m  <  n  <  (yc"11") 
or,  since  n  and  m  are  integers. 


(4.19) 


(4.20) 


(4.21) 

(4.22) 


(4.23) 


[  max]  .  m  +  1  <  n  <  [t  /cmin]  if  (T  /Cmax)  is  not  integral,  (4.24) 
and,  . 

[T  /Cmax]  -  m  <  n  <.  [T./C"110]  if  (T./Cwax)  is  integral. 

X  "  X  X 

It  should  be  noted  that  not  all  loop  equations  will  be  exactly  of 
the  above  form.  Different  pathologies  can  occur.  For  example,  the  loop  in 
Fig.  4.4  contains  links  in  opposing  directions.  Following  (4.15)  it  gives 
rise  to  this  equation: 


♦l  +  9i  +  ^  +  gh  +  ♦  *k  +  gk  +  lk  -  ♦  m+  9f  +  If  *  V 
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Essentially  as  we  traverse  a  loop  in  time  we  must  add  or  subtract  appropriate 
time  quantities  (according  to  the  direction  of  the  links  in  the  loop)  to 
return  to  the  starting  point  in  an  integral  number  of  cycle  times. 

4.3  Mixed-Integer  Linear  Programming  Formulation 

Lastly,  in  a  network  with  |N|  nodes  and  |L|  links,  exactly  |L|  -  |N|  +  1 
loop  constraint  equations  are  linearly  independent  [9],  A  corollary  to  this 
is  that  | N |  -  1  offset  variables  are  independent  and  they  must  correspond  to 
links  that  form  a  tree  subnetwork  of  the  general  network.  The  loops  corres¬ 
ponding  to  the  independent  set  of  constraints  are  called  a  fundamental  set  of 
loops  and  are  denoted  by 

The  objective  in  choosing  is  related  to  the  computational 

efficiency  of  mixed-integer  programming.  It  is  desirable  to  limit  the  number 
of  integer  lattice  points  of  the  problem,  i.e.,  maintain  close  upper  and 
lower  bounds  on  the  integer  variables.  This  is  achieved  by  choosing  loops 
that  contain  the  smallest  possible  number  of  links.  In  a  planar  network  this 
is  effected  by  choosing  the  meshes  of  the  planar  sketch  as  the  fundamental 
set  of  loops. 
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We  can  now  formulate  the  traffic  signal  network  coordination  problem 
as  the  following  mixed-integer  linear  programming  problem  (MILPP): 

MILPP:  Find  to 


min  7  f...  z.. 

(itj)  e  L  1J  " 

u 

subject  to:  z-  >  z^  (<j» - ^ ) 


l 


•ij  'Tij' 


k  =  1,2,...  ,ICjj 

V  (i.j)  e  L 


(i,j)  e  £ 


nU n*  i  < 

n£  integer. 


(♦ij  +  v  *  n»  c 


Vi  c  Z&f 


K-.  denotes  the  number  of  linear  segments  chosen  for  the  (i,j)  link  perform- 

w  1  U 

ance  function,  and  n£  and  n£  denote  respectively  the  lower  and  upper  bounds 
on  the  integer  variable. 


4.4  Test  Network  Solution 


The  following  data  describe  the  test  network  for  which  a  sample 
problem  is  actually  solved  on  the  IBM  370/165  using  IBM's  Mathematical 
Programming  Systems  (MPSX)  package.  Figure  4.5  is  a  sketch  of  the  test  net¬ 
work.  There  are  9  nodes  and  24  links  of  which  16  are  internal  to  the  net¬ 
work  and  8  are  input.  The  input  links  although  not  considered  in  this  prob¬ 
lem  will  be  mentioned  in  later  solutions.  Tables  4.1  and  4.2  present  the 
data  for  this  and  the  later  examples  of  the  test  network. 


In  this  example  only  offsets  are  decision  variables  with  green 
splits  and  cycle  time  fixed  according  to  Webster  [40].  The  critical  inter¬ 
section  is  No.  7  in  Fig,.  4.5,  yielding  an  80  second  cycle  time: 


_  1.5L  +  5  _  (1.5  x  9)  +  5  _ 
1-Y  ~  1-0.767 


79.5  sec  =  80  sec 


where  L  =  9  secs,  is  the  total  lost  time  at  that  intersection  (two  phases) 
and 


qi  q? 

+  y„  =  -r-  (link  123)  +  (link  ill)  =  0.767 
£  S1  s2 


Tables  4.3,  4.4  and  4.5  present  the  main  output  data  for  this  example. 
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Tdule  4.1 

Loops  of  the  Test  Network  and  Their  Corresponding  Integer  Variables 


Loop 

No. 

Integer 

Variable 

Link  Numbers 

1 

N1 

109 

110 

2 

N2 

103 

104 

N3 

105 

105 

4 

N4 

m 

112 

5 

N5 

101 

109 

104 

113 

6 

N6 

103 

111 

102 

114 

7 

N7 

115 

105 

113 

102 

8 

N8 

114 

106 

116 

108 

Table  4.2 


Link  Data  for  Test  Network* 


Link 

Travel 

Link 

Length 

Velocity 

Time 

Flow 

Sat.  Flow 

ilc . 

(ft) 

(ft/sec) 

(sec) 

(veh/sec) 

(veh/sec) 

101 

600 

32.27 

18.60 

.152 

.835 

102 

800 

32.27 

24.79 

.152 

.835 

103 

600 

32.27 

18.60 

.152 

.835 

104 

600 

41.07 

14.61 

.22 

.835 

105 

800 

32.27 

24.79 

.152 

.835 

106 

800 

41.07 

19.48 

.22 

.835 

107 

600 

41.07 

14.61 

.25 

.600 

108 

800 

41.07 

19.48 

.25 

.600 

109 

1000 

35.20 

28.41 

.175 

.500 

no 

1000 

41.07 

24.35 

.11 

.500 

111 

550 

35.20 

15.63 

.175 

.500 

112 

550 

41.07 

13.39 

.11 

.500 

113 

1000 

41.07 

24.35 

.12 

.360 

114 

550 

41.07 

13.39 

.12 

.360 

115 

1000 

35.20 

28.41 

.175 

.500 

116 

550 

35.20 

15.63 

.175 

.500 

Input 

Links 


117 

.175 

.500 

118 

.175 

.500 

119 

.152 

.835 

120 

.152 

.835 

121 

.12 

.360 

122 

.11 

.500 

123 

.25 

.600 

124 

.22 

.835 

(a)  Lost  time  for  each  link,  1  =  4.5 

secs. 

(b)  Platoon  length  on  input  links,  p 

=  1  cycle  (i.e. , 

flow 

to  be  continuous,  though  random  fluctuations  will  be  taken 
account  in  subsequent  formulations). 


PI atoon 

Length 

(cycle) 

.335 

.344 

.418 

.407 

.463 

.444 

.515 

.560 

.701 

.608 

.532 

.425 

.595 

.414 

.701 

.503 


assumed 

into 
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Table  4.3 


Test  Network  Results  and  Statistics* 

I .  Rows  1 61 

Columns  97 

Variables  258 

Integer  Variables  8 

Non-zero  Elements  709 

Density  1.70 


Time 

Iteration 

Node 

Functional 

(min) 

No. 

No. 

Value 

Continuous  Optimum 

.03 

7.6 

1 

43.3461 

First  Integer  Solution 

.06 

124 

18 

62.6695 

Optimal  Integer  Solution 

.07 

132 

21 

60.8741 

Optimality  Proved 

.21 

305 

57 

Time  of  Search 

.21 

305 

57 

III.  Number  of  Integer  Variables  Not  Integer  at  Continuous  Optimum  =  8 
Number  of  Integer  Solutions  Found  =  2 
Branches  Abandoned  While  Computing  =  50 


♦The  number  of  rows  and  columns  refer  to  the  linear  programming  coefficient 
matrix.  The  number  of  continuous  variables  includes  one  variable  for  each 
column  plus  one  slack  variable  for  each  row.  Thus  the  total  number  of 
variables  equals  the  number  of  rows  plus  columns.  Integer  variables  cor¬ 
respond  to  the  constraints  associated  with  each  network  loop.  The  density 
is  the  percentage  of  non-zero  elements  to  the  total  number  of  elements  in 
the  linear  programming  coefficient  matrix. 

In  solving  the  MILP,  the  integer  variables  are  first  treated  as 
continuous  and  the  resulting  LP  is  solved  yielding  the  continuous  optimum 
which  is  a  lower  bound  for  all  possible  feasible  integer  solutions.  Yet 
this  continuous  solution  is  an  infeasible  integer  solution.  At  each  level 
(Fig.  4.6)  an  integer  variable  is  chosen  with  a  non-integer  value  and  two 
problems  are  considered.  One  is  restricting  the  variable  to  be  greater  or 
equal  to  the  smallest  integer  greater  than  this  continuous  value,  the  other 
is  restricting  the  integer  variable  to  be  less  than  or  equal  to  the  largest 
integer  less  than  or  equal  to  the  continuous  value.  Then  each  of  these 
problems  is  solved  again  as  an  LP  and  the  next  node  considered  is  the  one 
that  produced  the  smallest  increase  in  the  objective  function.  This  process 
continues  until  we  finally  obtain  the  first  integer  solution.  The  program 
then  proceeds  to  follow  all  other  possible' paths  perhaps  finding  improved 
integer  solutions.  Once  all  feasible  possibilities  are  considered  optimality 
is  proved  and  the  last  integer  solution  is  chosen  as  the  optimal  solution  [7]. 
The  graphical  representation  of  this  process  is  called  a  Branch  and  Bound 
Diagram  (Fin.  4.6). 
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Table  4.4 


Integer  Nodes 


B-B  Node 

18 

21 

Functional 

Value 

62.6696 

60.8741 

Integer 

Variable 

Integer 

Integer 

N1 

1 

1 

N2 

1 

1 

N3 

1 

1 

N4 

1 

1 

N5 

-1 

-1 

N6 

-1 

-1 

N7 

0 

-1 

N8 

0 

0 

Figure  5.6  is  the  Branch  and  Bound  Diagram  corresponding  to  the 
MIL?  output  of  the  test  network.  The  nodes  are  numbered  according  to  the 
order  in  which  they  are  reached.  Nodes  with  x  indicate  fathoming*,  while 
proving  optimality.  Specific  nodes  and  values  were  omitted  from  the 
sketch  for  simplicity.  The  branching  is  indicated  along  the  line  segments 
connecting  two  nodes.  The  value  of  the  objective  function  is  given  beside 
each  node  while  an  asterisk  denotes  an  estimated  value. 


*(i.e.,  further  searching  along  these  branches  is  abandoned,  since  solution 
is  inferior  to  the  one  already  established). 


OPTIMAL  LP  'N  .... 
SOLUTION  )  43.35 


VN1  >  1 


45- 50  C  2 


3  ) 43.88 


N2  <  0/  \N2  >  1 


50.23(  s] 


4  ) 48. 69 


X  )51.12(  7 


N4  <  0, 


13  <  0 


6)55.57 


,N4  >  1 


stimated  value 


N5  <  -1 


X)  54.87(11 


N6  >  -1, 


15  >  0 


10) 60.75 


M6  <  -2 


53.22(12 


13)59.74  (X 


JI7  <  -1 


60.98  (  14 


[15  59.94  {  X 


NO  >  0/  Vl8  <  -1 

/  \  N5  >  -i 


kN5  <  -2 


62.56  (16 


N1  <  1 


62.67  (  18} 


(17)  (20J60.53  (19)61.58 
62.95  ^ 

V  J  \ 


50.77* 


feasible  62.66*  50787 
integer  optimal 

solution  integer 

solution 


Figure  4.6:  Branch  and  Bound  diagram. 

Splits  and  cycle  time  predeter¬ 
mined,  offsets  variable. 


Table  4.5 


Result  of  Computations* 

(C  =  80  sec;  cycle  time  and 
green  times  fixed  according  to  Webster) 


Link 

No. 

Offset 
sec  (cycle) 

Green  Times 
sec  (cycle) 

Delay 

SDF 

Degree  of 
Saturation 

101 

21.6  (.27) 
27.2  (.34) 

24.0  (.30) 

1.36 

0 

.607 

102 

24.8  (.31) 

1.36 

0 

.587 

103 

64.0  (.80) 

30.4  (.38) 

36.0 

0 

.479 

104 

16.0  (.20) 

31.2  (.39) 

.72 

0 

.676 

105 

47.2  (.59) 

31.2  (.39) 

11.28 

0 

.467 

106 

32.8  (.41) 

32.8  (.31) 

6.64 

0 

.850 

107 

16.8  (.21) 

39.2  (.49) 

1.04 

1.422 

.850 

108 

25.6  (.30) 

38.4  (.48) 

3.2 

1.574 

.868 

109 

26.4  (.33) 

40.8  (.51) 

10.88 

0 

.686 

110 

53.6  (.67) 

46.4  (.58) 

14.8 

0 

.603 

111 

32.8  (.41) 

32.0  (.40) 

9.6 

1.891 

.875 

112 

47.2  (.59) 

40.8  (.51) 

25.92 

0 

.686 

113 

16.8  (.21) 

45.6  (.57) 

5.36 

0 

.585 

114 

47.2  (.59) 

40.0  (.50) 

30.64 

.027 

.667 

115 

46.4  (.58) 

38.4  (.48) 

8.88 

.450 

.729 

116 

24.0  (.30) 

32.0  (.40) 

4.16 

1.891 

.875 

Input 

Links 


117 

46.4  (.58) 

12.16 

0 

.603 

118 

46.4  (.58) 

12.16 

0 

.603 

119 

24.0  (.30) 

24.56 

0 

.607 

1  rtA 

i  CV 

32.8  (.41) 

18.64 

0 

.444 

121 

31.2  (.39) 

22.16 

2.107 

.855 

122 

32,0  (.40) 

19.76 

0 

.550 

123 

38.4  (.48 

18.64 

1.574 

.868 

124 

30.4  (.38) 

21.52 

0 

.693 

♦The  delay  due  to  the  Saturation  Deterrence  Function  (SDF)  has  no  effect  on 
the  optimization  procedure  in  this  case  since  it  is  invariant  with  offsets, 
and  splits  and  cycle  time  are  kept  fixed.  The  same  is  true  for  the  delay 
values  due  to  the  LPF  on  the  input  links.  Both  factors  become  important 
when  splits  and  cycle  time  are  decision  variables  too  (see  Sections  5  and 
6). 


63 


5.  VARIABLE  SPLIT  FORMULATION 


In  this  chapter  we  add  the  split  at  each  node  of  the  network  as  a  decision 
variable  to  our  formulation  of  the  traffic  signal  network  optimization  pro¬ 
blem.  For  simplicity  we  consider  only  r.^  since  3  C-r^.  Now  the  objec¬ 
tive  function  of  delay  previously  dependent  on  offset  alone  is  a  function  of 
red  times  also.  Two  effects  are  added.  First,  we  must  consider  the  input 
links  of  the  network.  Second,  a  phenomenon  not  considered  before  has  an 
effect  on  the  split  decision— the  stochastic  behavior  of  traffic  flow.  This 
is  modeled  by  a  saturation  deterrence  function  (SDF).  Finally  we  give  the 
complete  formulation  followed  by  an  example. 


5.1  Link  Performance  Function  -  Approximations  by  Planes 

To  obtain  piecewise  linear  relations  in  the  (Y»r,z)  space  we  have 
to  fit  planes  by  selecting  triplets  of  points.  A  plane 


z  3  k,  +  k9  y  +  k,r 


(5.1) 


can  be  formed  from  the  points  Pj(y^ »r-,z.),  i 
relationships  from  geometry: 


1,2,3  by  the  following  simple 

\ 


triplets: 


Z1  - 

Vl  '  k3rl 

[<z2 

-  z,)(r3  -  r,)  - 

(z3  * 

■  z1><r2  - 

r^l/a 

C(*3 

"  Z-j )  (y2  “  Y-j  )  ~ 

(z2  * 

•  Z1)(t3  ' 

Y-j )]/° 

C(y2 

-  T,)(r3  -  r,)  - 

(y3  ' 

’  ,,l)(r2  ‘ 

r,)] 

g  to  Fig.  5.1  we  define  five  planes^ by  the 

K 

;  {P1’P2’P3} 

An 

:  "{P2 *P3P4 } 

Am 

:  {P.  |  z.=  0} 

A>I  V 

:  {P4*P5,P6  } 

:  ^P2*P4 *P6 } 

(5.2) 
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approximations  to  Link  Performance  Function. 


The  coordinates  are  determined  as  follows 


P1  :  Y1  =  "rmax 


rl  =  rmax 


Z1  =  2max(at  rl> 

P2  :  Y2  =  ’<P  '  W  =  C  *  rmax  "p 


r2  =  rmax 


Z,  =  Zm4„Ut  T.,) 
Z  min  2 


P  * 

r  o  • 


3  •  y3  -  y2 


r3  =  r4 


z3  =  z2 


p4  :  t4  =  0 


*A  *  C  -  P 


z4  =  ° 


P5  :  Y5  "  C  ~  rmin  'p 


r5  ~  rmin 


z5  =  0 
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P6  :  t6  =  p 
r6  =  r4 

z6  =  2nax(at  r4> 

rmax  1S  determined  by  capacity  requirements  as  follows: 

qp-gs  (5.3) 

or, 

OP  *  9m1n  5  *  (C  -  rmax>s- 

Consequently, 

rmax  *  C  '  w  (5’4) 

zm_„  at  a  certain  split  r  is  given  by  (3.27)  and  (3.36)  as  follows: 

§naX 

zmax  =  f(y  -  D  +  r  PO  -  y)  -  r 

r2  1 

zmax  =  pH  -  y>  >  r 

zmin  at  a  certain  split  r  is  given  by  (3.24): 

z  =  iB_-  J)f.  •  1  =  (P  ~  C  +  r)2  .  J 

min  2p  1  -  y  2p  1  -  y 

Planes  JLj  and  A,  i  are  parallel  to  the  r-axis,  according  to 
Theorem  3.5.  Therefore: 


and  we  can  determine  the  planes  by  two  points  P.(y*,z.)  only,  as  follows 
(see  Fig.  5.1):  1  1  1 


Al  :  |Pi  *^2} 

:  {p3,p4}. 
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The  plane  equations  in  this  case  are  of  the  form: 


z  =  k]  +  k2y 

Finally,  to  obtain  the  piecewise  linear  suface  as  a  function  of 
offsets  we  have  to  substitute  for  y  relationship  (2.2)  as  follows: 


Y  =  T  -  4» 

arrival  travel  offset 

time  time 

The  planesiij,... ,AV  enclose  a  convex  subspace  that  constitutes 

a  piecewise  linear  approximation  to  the  three-dimensional  surface  z(v,r). 

The  points  defining  the  planes  have  been  chosen  to  approximate  most  closely 
the  bottom  part  of  the  valley-like  enclosure,  i.e.,  where  the  optimal 
solutions  are  eventually  expected.  A  typical  set  of  z(y,r)  that  were  calcu¬ 
lated  for  a  rectangular-shapped  platoon  and  their  approximation  via  the 
planes  described  above,  are  shown  in  Fig.  5.2.  Note  that  the  line  connecting 
(P4,P2)  has  a  parabolic  shape,  according  to  equation  (3.24).  If  a  higher 

accuracy  is  warranted,  this  line  can  be  approximated  more  closely  by  a  series 
of  secant  planes,  rather  than  the  single  planefly. 


5.2  Adjusted  Traffic  Flow  Model 


In  Section  3  the  platoon  characteristics  were  calculated  according 
to  the  green  splits  at  the  intersection,  which  were  fixed  and  precalculated 
according  to  Webster's  formulas.  Since  in  the  current  context  the  upstream 
green  time  is  itself  a  decision  variable,  a  modified  procedure  has  to  be 
employed.  Given  the  cycle  time,  we  first  determine  nominal  splits  at  each 
node,  e.g.,  by  Webster's  method.  These  splits  are  used  to  Hotermine  the 
platoon  lengths  according  to  Section  3.2.  Though  the  final  splits  might  be 
slightly  different  than  those  assumed  for  calulating  p,  we  believe  that  this 
procedure  provides  a  satisfactory  approach.  An  alternative  approach  would  be 
to  conduct  an  iterative  process  in  which  platoon  legths  are  upaated  accord¬ 
ing  to  the  current  values  for  the  splits.  The  result  would  be  either  con¬ 
vergence,  or  a  "limit  cycle"  such  as  experienced  by  nonlinear  control  systems. 
In  the  first  case  no  distinguishable  difference  will  exist  between  split 
values  at  successive  iterations.  The  second  case  may  result  in  oscillations 
of  the  split  values  with  small  amplitudes  around  an  average  value.  Further 
computational  experimentation  is  required  to  evaluate  this  approach. 


5.3  Modeling  the  Network  Input  Links' 

In  the  current  formulation  splits  are  considered  as  decision  variables 
at  all  nodes  of  the  network.  Input  links  are  a  prime  determinant  of  the 
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fifes’  ffc 


ARRIVAL  TIME  (SEC) 


splits  at  the  signals  on  the  boundary  of  the  network  and  must  be  included  in 
the  program.  Traffic  on  these  links  is  assumed  to  be  independent  of  any 
neighboring  signals  and  therefore  will  be  modeled  as  a  continuous  flow,  i.e., 
with  p  =  1.  Though,  stochastic  effects  will  be  included,  as  explained 
below  in  Section  5.4. 


Figure  5.3:  Queue  evolvement  on  input  links  (p  =  r  +  g  =  C). 
Referring  to  Fig.  5.3,  for  p  =  r  +  g  =  C  we  have  for  the  clearance 

time  tpi 

qr  +  qtQ  -  stQ  =  0. 

Hence,  (5.5) 


The  total  delay  is  a  sum  of  two  areas: 

A  =  1  Qr2.  a  =  1  art  *  1 
”l  2  q*  *  m2  2  0  2  1  -  y 
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7  =  A  +  A  =  i  -9£ _ 

L  +  A2  2  1  -  y  ’ 


(5.6) 


and  tha  average  delay  is. 


Equation  (5.7)  is  equivalent  to  the  deterministic  delay  component  in  Webster's 
delay  formula[40]. 


5.4  Stochastic  Effects  -  The  Saturation  Deterrence  Function 

An  important  physical  phenomenon  is  missing  in  the  performance 
functions  of  most  signal  optimization  procedures.  At  flows  that  are  close  to 
capacity  but  still,  on  the  average,  below  it,  occasional  fluctuations  in  the 
size  of  the  platoon  can  lead  to  temporary  overflow  queues  that  seriously  de¬ 
grade  performance.  This  is  a  stochastic  phenomenon  which  has  a  consequence 
that,  as  average  flow  approaches  capacity,  average  delay  increases,  gradually 
at  first,  and  then  very  rapidly. 

A  representation  of  this  effect  is  needed  to  prevent  green  time 
from  approaching  its  lower  bound  too  closely.  It  is,  of  course,  possible 
to  put  a  sizeable  constraint  on  minimum  green  time  but  such  an  approach 
misses  the  main  idea  of  an  optimization,  which  is  to  tradeoff  delay  in  one 
part  of  the  system  against  that  in  another.  The  effect  is  particularly 
essential  if  cycle  length  is  a  variable  to  be  determined  by  the  optimization 
process  as  will  be  shown  later  in  Section  6.  This  is  because  of  the  tradeoff 
between  capacity  loss  at  short  cycles  and  the  inherently  large  delays  of 
long  cycles. 

Wormleighton  [43]  has  studied  the  effect  in  some  detail  motivated 
by  experience  with  the  Toronto  traffic  control  system.  Following  his  model, 
assume 

1)  Arriving  vehicles  come  in  platoons  or  other  periodic  function 
of  time  with  an  average  arrival  rate  at  time  t  of  q(t). 

2)  Arrivals  are  a  non-homogeneous  Poisson  process.  Thus  the 
number  of  vehicles  in  (t,t  +  C)  is  a  random  variable  having 
a  Poisson  distribution  with  mean 


t 


for  any  t. 
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3)  The  service  provided  by  the  green  time  is  deter¬ 
ministic  with  rate  s( vehicles/sec)  up  to  gs(vehicles/ 
cycle). 

If  the  state  of  the  intersection  is  examined  at  the  end  of  green, 
a  bulk  service  queuing  model  is  defined.  Let 

Q{0)  =  number  of  vehicles  in  queue  at  the  start  of  red  (end  of  green) 
—overflow  queue 

=  fC  =  average  number  of  vehicles  arriving  in  a  cycle 

S  =  gs  =  number  of  vehicles  that  can  be  served  in  a  green  time. 

x  =  Ac/s  =  fC/gs  =  utilization  factor  (degree  of  saturation). 

Wormleighton  finds  the  generating  function  of  Q(0)  and  calculates  E{Q(0) >  as 
a  function  of  S  and  x  over  S  e  [5,55]  and  x  e  [.2, .975].  Within  these  ranges 
E{Q(0) }  varies  from  essentially  zero  to  18  vehicles. 

We  wish  to  understand  how  the  presence  of  overflow  queues  affects 
delays  as  a  function  of  offset.  Insight  into  the  situation  is  obtained  by 
examining  the  deterministic  case.  A  theorem  is  proved  to  show  that  under 
certain  circumstances  delay  can  be  broken  into  two  components,  one  depending 
only  on  the  overflow  queue  and  the  other  being  the  delay  that  would  be  incur¬ 
red  in  the  absence  of  the  overflow  queue.  This  motivates  the  introduction  of 
the  overflow  queue  effect  into  the  link  performance  function  in  a  simple  way. 

Consider  the  sketch  of  queue  vs.  time  during  a  cycle  shown  in 

Fig.  5.4. 
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Let 


q(t)  =  arrival  rate  at  t (veh/sec) 

cc 

Ac  =  J  q(t)  dt  =  total  arrivals  in  cycle 
Q(t)  =  number  in  queue  at  time  t 
Q(0)  ~  overflow  queue  at  start  of  red 
Q0(t)  =  number  in  queue  at  t  if  overflow  queue 
had  been  Q{0)  *  0. 

Z  =  J  }(t)  dt  =  aggregate  delay  over  cycle  (veh/sec) 
0 

fc 

Zn  *  ,»  Qn(t)  dt  =  aggregate  delay  if  overflow 

0  **0  0 

queue  had  been  Q(0)  =  0(veh/sec) 

s  =  service  rate  during  green  while  vehicles  are 
present  in  queue (veh/sec) 


It  will  be  assumed  that  q(t)  <  s,  i.e.,  the  instantaneous  arrival  rate  never 
exceeds  the  instantaneous  service  rate. 

Theorem  5.1: 


If  Ac  =  9$,  i.e.,  \.he  number  of  arrivals  equals  the  maximum  number 
of  vehicles  that  can  be  served,  then 

(a)  the  overflow  queue  is  preserved  Q(C)  =  Q(0). 

(b)  the  aggregate  delay  is  composed  of  a  term  equal 
to  the  wait  of  the  overflow  queue  and  a  term 
equal  to  the  delay  in  the  absence  of  the  over¬ 
flow  queue: 
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(5.9) 


Z  =  Q(0)  C  +  ZQ 


Proof:  See  Fig.  5.4.  Because  AQ=  gs,  if  all  of  the  green  is  used  for  serving 
Vehicles,  the  queue  returns  to  its  starting  point  and  Q(C)  =  Q(0).  If  not  all 
the  green  is  used,  the  ending  queue  will  be  higher.  In  either  case 


Q(C)  >  Q(0) 


(5.10) 


Because  q(t)  is,  the  queue  never  increases  during  green  and  the  minimal  queue 
during  green  is  Q(C),  i.e.. 


Q(t)  £Q(C)  (t  in  green) 


(5.11) 


Therefore,  if  Q(C)  >  0,  some  queue  is  present  throughout  green  and  all  green 
time  is  used  for  serving.  In  this  case,  by  previous  argument,  the  queue  re¬ 
turns  to  its  starting  point 


Q(C)  =  Q(0) 


(5.12) 


On  the  other  hand  if  Q(C)  =  0,  Q(0)  *  0  by  (5.10),  and  (5.12)  is  trivially  true. 
Thus  in  any  case  (5.12)  holds,  proving  (a).  Notice  that  Q(t)  never  falls 
below  Q(0). 


Next  observe  that  the  order  of  service  does  not  affect  the  number 
in  queue  and  so  does  not  affect  the  aggregate  delays.  Suppose,  therefore, 
that  the  overflow  queue  vehicles  are  served  last.  Then  they  are,  in  fact, 
never  served  and  the  arriving  vehicles  are  served  just  as  if  Q(0)  =  0.  The 
delay  of  the  vehicles  arriving  during  the  cycle  is  Zq.  The  delay  of  the  over¬ 
flow  queue  vehicles  is  CQ(0).  The  total  delay  is  Z  =  Zn  +  CQ(0)  thereby  proving 
(b).  0 

We  use  the  above  result  to  argue  that  at  high  flow  whenAc^.  gs,  it 
is  reasonable  tc  separate  the  total  delay  into  two  components,  one  being  the 
result  of  the  presence  of  an  average  overflow  queue,  the  other  being  the  normal 
delay  of  vehicles  in  average  flow,  in  the  absence  of  an  overflow  queue. 

It  has  not  been  established  that  this  is  exactly  true  in  the  sto¬ 
chastic  case  and  so  the  term  involving  the  overflow  queue  should  be  regarded 
as  a  deterrence  function  which  involves  increasing  penalties  to  performance 
as  capacity  is  approached.  Thus  we  use  Wcrmleighton's  results  to  determine 
the  saturation  deterrence  function  (SDF),  which  enters  additively  into  the 
objective  function.  In  terms  of  delay  per  unit  of  time,  the  SDF  component  is 
simply  Q(0)  =  Q.  The  overflow  queue  Q  as  a  function  of  the  number  of  release 
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to 

in 


points  and  the  degree  of  saturation  is  listed  in  Table  5.1. 


5.5  Piecewise  Linearization  of  the  SDF 

We  must  be  able  to  represent  the  delay  of  the  SOF  in  a  form  amenable 
to  mixed  integer  linear  programming  techniques.  Although  there  is  no  simple 
analytic  expression  for  Q,  it  is  convex  in  r.  We  seek  ways  to  linearize  this 
convex  function.  Essentially  we  will  treat  this  delay  as  we  did  with  the  link 
performance  function,  first  using  secant  approximating  lines,  then  transfer¬ 
ring  them  into  the  constraint  set  with  corresponding  delays  contributing  addi- 
tively  to  the  original  objective  function  (  in  Section  4.3). 

Three  lines  will  be  used  to  estimate  this  delay  as  a  function  of  red 
times  (see  Fig.  5.5).  For  any  particular  line  we  restrict  r  to  be  such  that 
the  degree  of  saturation  remains  below  0.95. 


Table  5.2 


Example  of  Expected  Overflow  Queue  vs.  Red  Times 
C  =  60  sec 


s  =  .833  veh/sec 

S  -  gs 

f  “  .278  veh/sec 

r 

C 

£ 

C 

S 

X 

.05 

.95 

47.5 

.351 

0.0 

.10 

.90 

45.0 

.370 

0.0 

.15 

.85 

42.5 

.392 

0.0 

.20 

,80 

40.0 

.417 

0.0 

.25 

.75 

37.5 

.444 

0.0 

.30 

.70 

35.0 

.476 

0.0 

.35 

.65 

32.5 

.513 

0.0 

.40 

.60 

30.0 

.556 

0.0 

.45 

.55 

27.5 

.606 

.0i 

.50 

.50 

25. C 

.667 

.16 

.55 

.45 

22.5 

.741 

.38 

.60 

.40 

20.0 

.833 

1.2 

.65 

.35 

17.0 

,952 

7.48 

Piecewise  linearization  of 
Saturation  Deterrence  Function 


Our  capacity  constraint  reads 


Q  *  ♦  S  •  •  <  ^ •  »C 
a1J  u  1J 


or  in  terms  of  red  times 


r- .  £  C(1  -  f . ./s. -)* 
ij  v  ij'  ij' 


This  upper  limit  on  red  splits  is  obtained  by  considering 


°-95gijsij  1  V- 


Note  that  f. .C/g..s. .  is  the  degree  of  saturation  x.  The  other  two  lines  are 

'  J  '  J  '  J 

determined  by  the  following  procedure:  We  are  given  flows  and  saturation  flows 
for  every  link.  The  cycle  time  C  has  been  determined  according  to  Webster. 

For  the  first  line  we  choose  the  degree  of  saturation  x  r  .95  and  x  =  .90  which 
determines  the  corresponding  S  values  (S  =  gs  =  f  C/x,  i.e.,  the  number  of 
release  points).  Then  we  obtain  the  corresponding  Q  values  from  the  table 
either  directly  or  by  interpolation  if  necessary  for  each  link.  The  second 
line  is  determined  by  x  =  .85  and  x  =  .70,  yielding  again  directly  or  by  inter¬ 
polation  the  corresponding  Q  values.  An  example  to  the  final  approximation 
to  the  nonlinear  curve  is  given  in  Fig.  5.5.  The  data  for  the  original  curve 
is  given  ir.  Table  5.2. 

These  lines  appear  in  the  constraint  set.  We  require  the  $DF  delay, 
Q,  to  lie  above  each  line  and  represent  the  total  delay  from  this  effect  as 
the  sum  of  individual  delays  over  all  the  links  of  the  network.  It  should  be 
mentioned  that  the  first  line,  r..  <_  C(1  -  f../.95  s..)  is  not  handled  as  a 

IJ  •  J 

constraint  but  rather  as  an  upper  bound  to  the  red  split  variable.  Computa¬ 
tionally  this  is  more  efficient. 


5.6  Formulation  as  a  Mixed-Integer  Linear  Programming  Problem 

Now  we  are  able  to  present  the  complete  MILP  formulation  with  offsets 
and  splits  decision  variables.  We  denote  input  links  by  the  ordered  pair  (I,j) 
where  j  is  a  boundary  node  of  the  network,  to  distinguish  from  internal  links 
denoted  by  (i,j). 
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MILP :  Find  f..,  r..  to 

*  J  *  v 


min  I  +  Q.  ) 

(i,j)  e  L  1J 
(I.j)  e  L 


subject  to: 


zij  -2  (*ij’  r1j5 
zIj  —  zk(rt<) 


15  ‘ 


J'/ 


Qij  -  Q  (rij> 

QIj  -  ql(rIjJ 

(U)  e  +  V  =  n‘C 


k  =  1 ,K(LPF— internal  links) 
k'  =  1,...,K(LPF—  input  links) 

1*  =  1,...,L($DF— internal  links) 
1*  *  1,...,L(SDF— input  links) 

£ 


rij  "  ]j  =  rkj  "  ]j 


j  e  N;  i,k  e  P. 

J 


rij  -  (,j  +  aij)  =  9mj  +  1j  j  e  N;  i  c  Pjj  m  e  Pj 
rain 

r~  1  r-H  1  1  *  f^/0.95  s.,.  V  (i»j)  e  L 

(n£  integer) 


iJ  -  '  ij  -  '  ij'  1j 


"1  i  "ti  "t 


j  unrestricted  in  sign- 


5.7  Test  Network  Solution 

Our  test  network  of  Fig.  4.5  is  now  solved  for  offsets  and  splits 
variable.  The  cycle  time  is  again  determined  according  to  Wester  (C  =  80sec). 

The  results  are  given  in  Tables  5. 3,5.4,  and  5.5.  The  Branch  and  Bound  diagram 
corresponding  to  the  search  conducted  by  the  program  in  pursuit  of  the  optimal 
integer  values  is  given  in  Fig.  5.6.  (For  explanation  of  terms  see  Section  4.4). 


Table  5.3 


Test  Network  Results  and  Statistics 


1 .  Rows 

176 

Columns 

97 

Variables 

273 

Integer  Variables 

8 

Non- zero  Elements 

763 

Density 

1.58 

Time 

Iteration 

Node 

Functional 

(min) 

No. 

No. 

Value 

II.  Continuous  Optimum 

.04 

160 

1 

42.6181 

First  Integer  Solution 

.08 

236 

18 

61.4830 

Second  Integer  Solution 

.15 

250 

21 

59.6584 

Optimal  Integer  Solution 

.20 

384 

47 

59.6269 

Optimality  Proved 

.37 

576 

67 

Time  of  Search 

.37 

576 

67 

III.  Number  of  Integer  Variables  not  Integer  at  Continuous  Optimum  =  8 
Number  of  Integer  Solutions  Fpund  =  3 
Branches  Abandoned  While  Compute  rig  =  58 
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Figure  5.6  Branch  and  Bound  diagram.  Cycle  time  predetermined,  offsets 
and  splits  variable. 


Table  5.4 
Integer  Nodes 


B-B  Node 

18 

21 

27 

Functional 

Value 

61.4830 

59.6584 

59.6269 

Integer 

Variable 

Integer 

Integer 

Inteqer 

N1 

1 

1 

0 

N2 

1 

1 

1 

N3 

1 

1 

1 

N4 

1 

1 

1 

N5 

-I 

-1 

1 

N6 

-1 

-1 

-1 

N7 

0 

-1 

-1 

N8 

0 

0 

0 

.(The  B-B  nodes  refer  to  Fig.  5.6). 
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Table  5.5 


| 

I 


Results  of  Computations 
(C  *  SO  sec.  according  to  Webster) 


1 

I 

'i 

Link 

No. 

Offset 

sec  (cycle) 

Green  Times 

sec  (cycle) 

Delay 

Delay 
Due 
to  SDF 

Degree  of 
Saturation 

if 

101 

24.0  (.30) 

21.6  (.27) 

2.72 

0 

.674 

! 

102 

28.8  (.36) 

23.2  (.29) 

2.08 

0 

.628 

’• 

103 

62.4  (.78) 

29.6  (.37) 

34.64 

0 

.492 

f 

104 

17.6  (.22) 

29.6  (.37) 

1.36 

0 

.712 

105 

46.4  (.58) 

29.6  (.37) 

10.64 

0 

.492 

106 

33.6  (.42) 

30.4  (.38) 

7.2 

0 

.693 

k 

107 

14.4  (.18) 

40.8  (.51) 

.08 

1.067 

.817 

Uf 

108 

24.8  (.31) 

40.0  (.50) 

2.56 

1.339 

.833 

£ 

109 

18.4  (.23) 

39.2  (.49) 

13.36 

.279 

.714 

h 

110 

-18.4  (-.23) 

49.6  (.62) 

13.52 

0 

.355 

P 

111 

35.2  (.44) 

33.6  (.42) 

11.92 

1.505 

.833 

1 

112 

44.8  (.56) 

39.2  (.49) 

23.6 

0 

.449 

I 

113 

24.0  (.30) 

47.2  (.59) 

0.0 

0 

.565 

ft 

f 

114 

44.0  (.55) 

49.8  (.51) 

27.36 

0 

.654 

115 

41.6  (.52) 

40.0  (.50) 

8.32 

0 

.700 

| 

116 

24.8  (.31) 

31.2  (.39) 

4.48 

2.035 

.897 

Input 

Links 


117 

49.6  (.62) 

9.92 

0 

.565 

118 

49.6  (.62) 

9.92 

0 

.565 

119 

21.6  (.27) 

26.4 

0 

.674 

120 

30.4  (.38) 

20.08 

0 

.479 

121 

29.6  (.37) 

23.04 

2.410 

.901 

122 

33.6  (.42) 

18.64 

0 

.524 

123 

37.6  (.47) 

19.36 

1.784 

.887 

124 

29.6  (.37) 

21.92 

0 

.712 
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Section  5.2  indicates  that  a  cause-effect  relationship  exists  among 
the  signal  splits  which  are  decision  variables  of  the  program,  and  the  platoon 
model  of  traffic  flow,  which  is  an  input  parameter  to  the  optimization  pro¬ 
cedure.  Iterating  through  the  sequence:  platoon  -*  splits —platoon  — spl its, 
the  optimized  splits  were  used  to  recalculate  platoons.  The  results  are  given 
in  Table  5.6.  The  initially  calculated  splits  are  satisfactory,  in  the  sense 
that  no  appreciable  differences  in  the  resulting  platoons  are  produced  by 
iterating. 


Table  5.6 


Iteration  on  Platoon  Length 


(C  =  80  sec) 


Link  No. 

»1 

P2 

A 

101 

.335 

.27 

.312 

.27 

102 

.344 

.29 

.279 

.28 

103 

.418 

.37 

.398 

.37 

io4 

.407 

.37 

.398 

.37 

105 

.463 

.37 

.434 

.37 

106 

.444 

.38 

.423 

.37 

107 

.515 

.51 

.497 

.47 

108 

.560 

.50 

.582 

.47 

109 

.701 

.49 

.745 

■  .49 

n 

10 

.608 

.62 

.588 

.62 

11 

!1 

.532 

•  42 

.515  ’ 

.42 

ll 

12 

,425 

.49 

.442 

.49 

ii 

13 

.695 

.59 

.618 

*61 

11 

14 

•414 

.51 

.393 

.51 

n 

15 

*70 

.§  0 

,746 

.51 

116 

.503 

.39 

.530 

.41 
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This  can  be  easily  interpreted  with  reference  to  Fig.  5.7.  Consider¬ 
ing  for  example  a  two-phase  intersection  we  have: 


g-j  +  g£  +  21  =  C 

split  g  varies  from  a  minimum  of  gmin  (determined  by  capacity  requirements)  and 
gmax  (determined  by  capacity  requirements  in  the  conflicting  phase).  Therefore, 

max  _  p  ai  _  min 
g1  -  C  -  21  -  g2 


The  combination  of  the  SDF  values  for  both  approaches  will  have  the  parabolic 
shape  shown  in  Fig.  5.7.  Consequently  any  choice  of  offsets  will  most  likely 
lie  in  the  neighborhood  of  the  minmum,  which  is  very  closely  approximated  by 
Webster's  split-decision  formula. 


Overflow  Queue 


0  gmin  gmax 

Fig.  5.7:  Effect  on  overflow  queue  of  split  value. 
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6.  The  Network  Synchronization  Problem 


In  the  preceding  sections  we  have  formulated  the  traffic  signal  network 
problem  as  minimizing  delay  due  to  the  deterministic  link  performance  function 
(LPF)  and  the  stochastic  saturation  deterrence  function  (SDF)  subject  to  loop 
and  other  physical  constraints.  We  first  considered  the  offset  optimization 
problem  (the  network  coordination  problem)  and  then  added  the  splits  at  each 
intersection  as  decision  variables.  How  we  shall  let  the  cycle  time,  common 
for  all  the  intersections  of  the  network,  be  a  decision  variable  too.  This 
will  be  termed  as  the  network  synchronization  problem. 


6.1  The  Optimal  Cycle  Time  in  a  Network 

Experience  of  researchers  and  practitioners  in  the  urban  traffic  con¬ 
trol  field  has  shown  that  the  cycle  time  may  well  be  the  most  important  among 
the  decision  variable  [3,17,38].  Not  only  does  it  tie  together  the  offsets  in 
conjuction  with  integer  variables  but  it  also  provides  for  the  necessary  capa¬ 
city  to  serve  the  traffic  demand  at  each  intersection.  The  net  serving  capa¬ 
bility  of  an  intersection  is  1  inearly  proportional  to  the  sum  of  the  effective 
green  times  over  all  the  phases  of  the  intersection,  i.e.,  capacity  is  propor¬ 
tional  to 


(6.1) 


Since  L,  the  total  lost  time,  is  a  fixed  quantity  at  a  particular  intersection, 
the  net  capacity  increases  with  cycle  time.  For  example,  for  a  two  phase  inter¬ 
section  with  a  9  second  total  lost  time  (4.5  seconds  per  phase),  capacity  in¬ 
creases  by  varying  cycle  time  from  40  seconds  to  120  seconds  by  the  difference 
between  9/40  and  9/120,  i.e.  by  15%.  This  meager  difference  in  net  capacity 
may  have  a  decisive  impact  on  delay,  particularly  at  high  degrees  of  saturation. 
Similar  characteristics  are  well  known  in  many  other  queueing  processes  [  29  ]. 
On  the  one  hand  an  increase  in  cycle  time  is  usually  followed  by  an  increase 
in  red  times  on  individual  phases  and  consequently  by  an  increase  in  the  aver¬ 
age  waiting  time  as  expressed  by  the  LPF.  On  the  other  hand,  a  decrease  in 
cycle  time  reduces  capacity,  thus  precipitating  additional  delay  due  to  the 
randomness  in  arrivals  of  vehicles.  This  effect  may  become  particularly  severe 
when  flow  approaches  saturation  and  is  modeled  by  the  SDF. 

The  interplay  between  the  LPF  and  the  SDF  in  the  objective  function 
for  variable  cycle  time  was  studied  with  respect  to  the  test  network  intro¬ 
duced  in  Section  4.4.  The  MILPP  formulation  of  Section  5.7  was  run  for  various 
f i xed  cycle  times  in  the  range  of  40  seconds  to  120  seconds  to  optimize  offsets 
and  splits.  The  resultant  values  *or  the  objective  function,  broker  down  into 
its  two  components,  are  given  in  Table  6.1.  A  graphical  illustration  of  the 
relationships  is  presented  in  Fig.  6.1.  It  becomes  evident  that  the  optimal 
cycle  time  for  the  network  constitutes  a  least-cost  equilibrium  point  between 
delays  caused  by  deterministic  effects  and  delays  contributed  by  stochastic 
effects.  While  the  first  delays  usually  increase  with  cycle  length  (though  at 
a  decreasing  rate),  the  latter  decrease  with  it  owing  to  the  decrease  in  the 
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Table  6.1 


Variation  of  Delay  with  Cycle  Time 
for  Test  Network  of  Section  4.4[veh  x  sec/sec] 

CYCLE  TOTAL 

(sec)  LPF  SDF  DELAY 


50 

27.5 

55 

31.3 

63.7 

37.7 

70 

42.5 

80 

49.2 

90 

53.5 

103 

57.1 

120 

59.8 

43.6 

71.1 

29.5 

60.8 

15.4 

53.1 

13.1 

55.6 

10.4 

59.6 

9.4 

62.9 

8.7 

65.8 

8.4 

68.2 

WEBSTER 


degree  of  saturation  (the  load  factor).  They  approach  asymptotically  from  above 
the  minimal  cycle  time  for  the  network,  which  is  the  theoretical  minimal  cycle 
time  for  the  most  heavily  loaded  intersection  if  all  flows  were  strictly  deter¬ 
ministic.  These  characteristics  are  in  complete  analogy  with  the  behaviour  of 
delay  with  respect  to  cycle  time  at  a  single  intersection  and  were  studied  by 
Wardrop[39  ],  Webster[40  ]  and  other  researchers  [26,32].  However,  the  impli¬ 
cations  regarding  signal  settings  in  a  network  are  different  and  a  single 
intersection  analysis  would  virtually  never  give  the  optimum  settings  for  the 
network. 


6.2  Variable  Cycle  Formulation 

As  shown  above,  to  optimize  signal  settings  in  a  network,  any  optimi¬ 
zation  procedure  must  make  a  judicious  choice  of  cycle  time.  To  introduce  the 
cycle  time  as  a  decision  variable  our  formulation  of  the  network 
problem  has  to  be  modified.  It  will  be  observed  that  by  a  change  of  variable 
the  cycle  time  enters  linearly  into  the  LPF  and  only  with  slight  modifications 
into  the  SDF.  Examination  of  the  loop  constraints  with  C  variable  will  provide 
the  motivation  for  this  change  of  variable.  The  input  links  will  not  be 
affected.  Special  consideration  must  be  given  to  the  transformed  objective 
function  and  its  linearized  representation. 


For  variable  cycle  time  our  objectives  become  twofold.  First  we  must 
examine  our  model  of  delay  and  the  physical  constraints  with  this  added  variable. 
Then  we  must  insure  that  the  optimization  procedure  can  handle  these  changes. 
Examining  the  loop  constraints (2. 5)  we  observe  that 

(Id)  «  t  *«  *  «1C 


for  variable  cycle  time  there  is  a  nonlinearity  in  the  term  n£C.  Common  means 
of  handling  such  a  nonlinearity  as  separation  of  variables  or  Taylor  approxi¬ 
mation  do  not  seem  to  be  of  use  in  this  case  since  nfc  is  an  integer  variable. 

Yet  by  dividing  each  term  of  this  equation  by  C  we  can  remove  this  nonlinearity, 
introducing  terms  like  $/C,  r/C  and  1/C  which  are  also  nonlinear.  This  suggests 
a  change  of  variable  as  follows 


=  <J>/C,  r'  =  r/C,  w  =  1/C 


(6.2) 


i.e.,  we  transform  variables  from  seconds  to  dimensionless  fractions  of  a  cycle 
time.  From  this  point  on  we  drop  all  primes  on  variables,  understanding  that 
they  are  measured  in  fractions  of  a  cycle  time. 

So  far  platoon  lengths  p  have  been  determined  in  seconds.  Following 
the  above  transformation,  p  will  be  measured  in  fractions  of  a  c)Je  too.  We 
assume  that  the  ratio  of  platoon  length  to  cycle  length  (in  seconds)  remains 
invariant  as  we  allow  the  cycle  time  to  vary.  This  assumption  is  reasonable 
since  the  platoon  length  on  any  link  (i,j)  is  considered  to  be  proportional  to 
the  upstream  green  which  varies  linearly  with  cycle  time.  Thus  the  physical 
constraints  remain  linear  wiht  the  above  transformation.  Now  we  must  examine 
the  constraints  corresponding  to  the  different  components  of  the  delay  function. 
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The  link  performance  function  has  been  expressed  in  terms  of  arrival 
times  y  and  splits,  according  to  (2.2) 

arrival  time  offset  travel  time 

y  *  ♦  -  t 

Transforming  arrival  time  to  fractions  of  a  cycle  we  obtain 

Y *  =  <#>  *  -  tw  (6.3) 

so  the  reciprocal  cycle  time  enters  linearly  into  the  constraints  representing 
the  LPF.  The  LPF  for  the  input  links  is  independent  of  C  and  requires  no  such 
transformations. 

6.3  SDF  -  Approximation  by  Planes 

Diecycle  time  will  now  be  introduced  intotheSDF  through  its  reciprocal 
w  *  i/C.  Referring  to  Table  5.1  we  notice  that  the  expected  overflow  queue  Q 
in  a  function  of  the  degree  of  saturation  x  and  the  number  of  release  points 
in  the  cycle,  S  (i.e.,  the  capacity  of  the  approach).  By  definition  we  have. 


s  =  sgC  = 

(6.4) 

(g  and  r  are  given  in  fractions  of  cycle  time) 

..  _  f  _  f 
x  sg  s(l  -  r) 

(6.5) 

Solving  for  r  and  w  we  obtain 

(6.6) 

II 

X 

(6.7) 

We  may  now  transform  the  triplet  of  points  (Q,q,x)  into  (Q,r,w)  which  -coincides 
with  our  decision  space. 

Figure  6.2  shows  a  typical  relationship  between  Q,r,  and  w  which  sug¬ 
gests  planar  approximations  as  before.  The  first  plane,  an  upper  bound  on  red 
times,  confining  the  degree  of  saturation  to  values  below  0.95,  will  be 


r  <  1 


f 

0.95s 


(6.8). 


as  in  the  case  with  fixed  cycle  time  (Section  5.5).  The  other  two  planes  will 
be  determined  by  the  triplets  of  points: 
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Variation  of  overflow  queue 
with  split  and  cycle  time. 


W=  V60  sec1 
W  =  Vfoo  sec1- 


Plane  I: 


Plane  II: 


a 

S 

X 

7.61 

13 

0.95 

2.81 

i5 

0.90 

2.11 

35 

0.90 

1.76 

15 

0.85 

1.23 

35 

0.85 

0.37 

15 

0.70 

The  values  for  Plane  II  were  obtained  by  interpolation  from  the  values  given 
in  Table  5.1.  Note  that  the  values  for  Q  change  more  pronouncedly  on  the 
horizontal  axis  in  Table  5.1  than  in  the  vertical  axis.  This  means  that  Q  is 
more  sensitive  to  variations  in  split  than  in  cycle  time,  as  evidenced  also  by 
the  plot  in  Fig.  6.2.  As  a  result  the  planes  representing  the  SDF  (appearing 
as  contraints  in  the  MILP  formulation) 


Q  +  er  +  bw  >0 


(6.S) 


show  that  the  coefficients  a  are  ccnslderal  iy  larger  in  magnitude  than  the 
coefficients  b. 

6.4  Revised  Objective  Function 

Introducing  C  as  a  variable  has  the  most  significant  effect  on  the 
objective  function.  Prior  to  this  change  -cur  objective  function  was  the  simple 
linear  expression 


( 


I  +  Q«h  where*,. 

i,j)  1J  1J 


is  in  seconds.  Incorporating  C  as  a  decision  variab<f,  thus  transforming  vari¬ 
ables  to  dimensionless  .quantities,  introduces  a  chance  ir<  the  objective  function. 
Now  we  must  consider  z-j  in  fractions  of  a  cycle.  We  define  z*  "  z/C,  ov 


z  =  Cz*  =  £ 


(6.10) 


Thus  our  objective  function  now  reads 

mi"  I  (t'-.— 51+ 

(i.j)  w 

The  transformation  of  z  to  dimensionless  z',  provides  us  with  a  LPF  that  is 
Invariant  with  cycle  time.  No  change  occurs  in  the  variable  Q..  due  to  the 
SDF. 
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We  cannot  proceed  in  a  straightforward  manner  since  z/w  is  nonlinear 
(we  drop  agiin  primes).  Although  it  can  be  shown  that  z/w  is  quasi-convex, 
which  may  lend  itself  to  gradient  techniques,  it  is  not  convex.  In  order  to 
maintain  the  MILP  formulation,  we  introduce  again  piecewise  planar  approxima¬ 
tions  of  z/w.  That  is,  we  approximate 


(6.11). 


by  a  set  of  plane*  forming  a  convex  region.  As  before,  the  planes  will  be 
transferred  into  the  constraint  set  of  the  objective  function,  which  will  now 
read 


min 


I 

(i»j) 


(d 


1j 


+  Q. .). 
13 


Care  must  be  taken  in  approximating  a  non-convex  surface  by  a  set  of 
planes  forming  a  convex  surface.  The  difficuHy  is  that  we  may  cut  away  large 
portions  of  the  feasible  region.  Thus  we  want  to  approximate  the  surface  as 
closely  as  possible  in  the  domain  in  which  we  expect  the  solution  to  lie  and 
only  exclude  extreme  portions  of  the  feasible  region  where  we  are  quite  sure 
will  be  unlikely  for  a  solution  to  occur.  This  is  exactly  the  case  in  the 
approximation  to  the  LPF,  wherp  there  is  a  non-convex  portion  corresponding  to 
very  high  green  splits. 

The  following  three  planes  are  chosen  to  approximate  the  d  =  z/w  sur¬ 
face  for  internal  links  of  the  network.  Each  of  the  planes  is  determined  by 
three  points  P  (see  Fig.  6.3). 


z 

w 

d  «  z/w 

h 

* 

• 

0.15 

Q.018 

0.15/0.018 

Plane  I: 

P2 

* 

* 

0.15 

0.0.4 

0.15/0.014 

p3 

• 

• 

0.00 

1/120 

0.00 

P1 

• 

• 

0.15 

0.01 

0.15/0.01 

Plana  II: 

P2 

• 

• 

0.15 

0.014 

0.15/0.014 

P3 

• 

* 

0.00 

1/120 

0.00 

P1 

• 

• 

0.16  ' 

0.018 

0.15/0.018 

Plane  III: 

P2 

* 

• 

0.00 

0.00 

0.00 

P-. 

• 

0.00 

0,00 

0.00 
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Figure  6,3:  Linear  approximations 


The  planes  have  been  constructed  taking  into  account  that  the  lower  limit  on 
w  is  1/120  (corresponding  to  an  upper  limit  on  cycle  time  of  120  seconds), 
z  =  0.15  has  been  chosen  as  the  pivot  value,  representing  a  typical  average 
delay  per  car  figure  (in  fractions  of  cycle  time). 

The  planes  approximating. the  LPF  on  input  links  to  the  network  have 
been  constructed  similarly,  with  a  pivot  value  of  z  =  0.25.  The  reason  is  that 
platoons  on  these  links  are  assumed  to  have  full  cycle  length,  i.e.,  flow  is 
uniformly  distributed  over  time.  Thus  the  average  delay  for  vehicles  on  those 
links  is  expected  to  be  higher  than  on  the  internal  links  of  the  network. 

S.5  Mixed-Integer  Linear  Programming  Formulation 

The  complete  optimization  problem  including  simultaneously  offsets, 
splits  and  frequency  (cycle  time)  as  decision  variables  is  formulated  as  follows: 


MILPP:  Find  r^r  w  to 

min  X  ^ii^ii  +  Oh) 
(i,j)  cl  1J  1J 


subject  to: 


Z1j^ 

i!t<W *> 

k  «  1,. 

.-.  ,K(LPF) 

Qij  - 

Q  '(rH.w) 

1  *  1.. 

. .,L(SDF) 

dij  £ 

**u-> 

m  =  1 ,. 

..^(representation 

X 

(i  J) 

c  .<♦«  +  V  *  "* 

i 

rij  ’ 

V  =  rkj  '  V 

j  e  N; 

1,k  «  Pj 

rij  " 

(’id  +  aij>"  =  Vj  + 

S  e  N; 

i  e  Pji  j  e  Pj 

jnin 
r  ' 

S  rfj  <  1  -  f^/0.95  ,y 

(i »0 )  £ 

:  L 

/1n  <  »</** 

nj  <  n,  £  nj  (n(  Integer) 

unrestricted  in  sign. 


6.6  Test  Network  Solution 


The  test  network  of  Fig.' 4.5  is  now  solved  for  offsets,  splits  and  cycle 
time  variable  simultaneously.  The  results  are  given  in  Tables  6. 2,6. 3,  and  6.4 
below.  The  Branch  and  Bound  diagram  followed  by  the  MILP  program  in  search  for 
the  optimal  integer  set  is  illustrated  in  Fig.  6.4.  As  seen  in  Table  6.3  the 
optimal  solution  in  this  case  is  a  considerable  improvement  over  those  obtained 
previously  (see  Section  4.4  and  Section  5.7),  due  to  the  decisive  role  played 
by  the  cycle  time  in  the  total  optimization.  By  considering  all  the  decision 
variables  simultaneously  (objective  function  =  53.80),  network  performance  was 
improved  by  11.6%  with  respect  to  settings  obtained  when  the  decision  variables 
were  considered  sequentially  (objective  function  *  60.87). 


!  Table  6.2 

j  - 

! 

■  Test  Network  Results  and  Statistics 


I .  Rows 

248 

Columns 

121 

Variables 

369 

Integer  Variables 

8 

Elements 

1027 

Density 

1.12 

*  Time 
(min) 


Iteration 

No. 


Node  Functional 
No.  Value 


Continuous  Optimum 

.10 

230 

First  Integer  Solution 

.24 

376 

Optimal  Integer  Solution 

.25 

385 

Optimality  Proved 

.40 

499 

Tima  of  Search 

.40 

499 

1 

40.0006 

16 

59.6186 

17 

53.8010 

23 

23 

III.  Number  of  Integer  Variables  Not  Integer  at  Continuous  Optimum  *  8 
Number  of  Integer  Solutions  Found  *  2 
Branches  Abandoned  While  Computing  *  12 
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Table  6.3 


Integer  Nodes  and  Functional  Values 


B-B  Node 

16 

17 

Functional 

Value 

59.6186 

53.8010 

Frequency 

wCsec"1) 

0.01565 

0.01569 

Integer 

Variable 

Integer 

Integer 

N1 

1 

1 

N2 

1 

1 

N3 

1 

1 

N4 

1 

1 

N5 

-1 

-1 

N6 

-1 

•  -1 

N7 

0  . 

■0 

N8 

-1 

0 

(The  B-B  Nodes  refer  to  Fig.  6.4) 

In  order  to  check  the  quality  of  our  planar  approximations  to  the 
objective  function  d  =  z/w,  the  test  network  was  re-solved  with  the  optimal 

value  of  w  =  0.01569  sec“^(C  =  63.8  sec),  kept  fixed.  The  nonlinearity  of  the 
objective  function  is  circumvented  in  this  case.  The  same  integer  set  resulted, 
with  all  other  variables  being  very  close  to  the  variable-cycle  optimal  solution. 
The  functional  value  value  was  53.1272,  i.e.,  a  difference  of  only  U25%. 

6.7  Sensitivity  Analysis 

Since  in  practice  flows  change  frequently,  the  test  network  was 
further  used  to  analyse  the  sensitivity  of  delay  with  respect  to  flows.  All 
flows  in  the  network  were  changed  by  the  same  percentage  and  the  MILPP  formu¬ 
lation  of  Section  5.7  was  used  to  optimize  offsets  and  splits  for  various  cycle 
times.  The  results  are  shown  in  Fig.  6.5.  The  decisive  role  played  by  the 
cycle  time  in  reaching  optimum  operating  conditions  is  illustrated  here  even 
more  emphatically  than  before.  Again,  similarity  with  single  intersection 
behavior  is  apparent[26,40]. Aclear  conclusion  is  that  it  certainly  makes  good 
sense  to  try  and  adapt  cycle  length  to  actual  traffic  conditions.  However,  how 
to  do  this  in  an  optimal  manner,  taking  into  account  the  transient  phenomena 
involved  in  this  process,  requires  further  research. 


Table  6.4 


Results  of  Computations 


Link 

No. 

Offset 

sec  (cycle) 

Green  Times 

sec  (cycle) 

Delay 

(LPF) 

Delay 
Due  to 
SDF 

101 

28.0  (.44) 

17.2  (.27) 

4.84 

0 

102 

29.3  (.46) 

17.2  (.27) 

2.23 

0 

103 

38.2 

1.60] 

24.2  (.38) 

12.43 

0 

104 

25.5 

;.40] 

21.7  (.34 

5.99 

.943 

105 

32.5 

[•51  ] 

21.7  (.34) 

3.95 

0 

106 

31.2 

[.49] 

24.2  (.38) 

5.74 

0 

107 

21.0 

[.33] 

30.6  (.48) 

3.12 

1.660 

108 

25.5 

[-40] 

30.0  (.47) 

3.00 

1.930 

109 

44.6 

[-70] 

28.7  (.45)  ' 

8.09 

1.255 

110 

19.1 

[.30] 

37.6  (.59) 

2.87 

0 

111 

.  33.8 

[.53] 

24.9  (.39) 

12.43 

2.244 

112 

30.0 

(.47] 

28.7  (.45) 

4.94 

0 

113 

33.8 

[.53, 

37.6  (.59) 

4.72 

0 

114 

30.6 

(.48 

33.1  (.52) 

14.60 

0 

115 

24.9 

(.39, 

30.6  (.48 

10.26 

.690 

116 

24.2 

(.38 

24.9  (.39) 

4.14 

2.244 

Input 

Links 


117 

37.6  (.59] 

9.43 

0 

118 

33.1  (.52) 

12.11 

0 

119 

21.7  (.34) 

18.04 

0 

120 

24.2  (.38) 

16.12 

0 

121 

23.6  (.37) 

18.29 

2.714 

122 

24.9  (.39) 

16.06 

0 

123 

30.0  (.47) 

15.42 

1.930 

124 

24.2  (.38) 

17.14 

0 

Degree  of 
Saturation 

.674 

/674 

.479 

.775 

.535 

.693 

.868 

-.887 

.778 

.373 

.897 

.489 

.565 

.641 

.729 

.897 


.593 

.673 

.535 

.479 

.901 

.564 

.887 

.693 
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Branch  and  Bound  diagram.  Offsets,  split  and  cycle 
time  variable. 


respect  to  flows. 


7.  Applications  and  Computational  Results  of  MITROP 

This  chapter  presents  computational  results  that  were  obtained  with  the 
mixed-integer  traffic  optimization  program  (MITROP)  for  the  UTCS-1  Test  Net¬ 
work,  and  for  an  arterial  in  Waltham,  Massachusetts. 

These  results  are  considered  preliminary  since  it  is  the  belief  of  the 
authors  that  further  iterations  with  improved  data  are  necessary  to  obtain 
results  that  will,  in  fact,  optimize  the  performance  of  the  actual  networks. 

In  contrast  with  previous  sequential  optimization  methods,  MITROP  applies  a 
rigorous  optimization  procedure  for  determining  all  the  decision  variables  in 
a  signal-controlled  network— off sets ,  splits,  and  cycle  time— simultaneously. 
It  is  of  great  importance  that  the  data  truly  reflect  the  physical  situation 
to  which  it  is  applied.  In  particular  the  most  heavily  loaded  intersections 
should  be  examined  carefully  in  order  to  obtain  accurate  results  [35]. 


7.1  UTCS-1  Test  Network 

Here  we  present  some  computational  results  for  the  UTCS-1  Test  Net¬ 
work  (Fig.  7.1).  Every  traffic  network  has  its  own  peculiarities  in  topology 
and  certain  care  must  be  taken  in  preparing  the  input  data  before  an  optimi¬ 
zation  procedure  is  employed.  We  shall  indicate  how  data  is  transformed  and 
in  particular  the  handling  of  extra  phases  as  auxiliary  links.' 

A.  Flows 

For  any  given  link  we  are  given  primary  flows  secondary  flows,  and 

turning  movements.  First,  we  add  primary  and  secondary  flows  (f  +  f  )  .to 

P  ^ 

obtain  a  total  flow  fT.  This  must  then  be  adjusted  for  turning  movements. 

Given  percentages  of  turning  movements,  left(%).=  a,  thru  (%)  =  b,  right 
(%)  =  c  such  that  a  +  b  +  c  =  1  we  adjust  as  follows: 


fT(l  +  3/4  a  +  l/4c)  =  F  (7.1) 


This  is  again  adjusted  for  trucks  and  buses,  i.e.,  3%  increase  for. buses 
(except  on  K-street  where  we  use  6%)  and  a  5%  increment  for  trucks^.  So  the 
final  total  flow  is  1.08F  or  1.11F  on  K-street. 

Platoon  lengths  are  computed  as  described  in -Section  3.2. 

B.  Extra  Phases 


Additional  physical  node  constraints  may  occur  in  a  network  due  to 
extra  phases.  For  instance,  or.  ai.y  link  an  extra  phase,  i.e.,  an  additional 
phase  to  accommodate  turning  movement,  may  be  modeled  as  an  auxiliary  link 
running  parallel  *;o  the  main  link.  This  may  introduce  two  additional  types 
of  constraints,  one  reflecting  capacity  and  the  other  sequencing. 
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If  on  link  1  we  have  a  phase  for  straight  thru  traffic  plus  a  phase 
for  left  turning  movements  which  partially  coincides  with  the  main  phase  in 
time,  an  extra  link  j  incident  to  node  n  and  In  the  direction  of  link  1  may  be 
Introduced. (Fig.  7.2). 


Figure  7.2  Schematic  of  an  extra  phase  j. 


For  the  duration  of  green  in  direction  i,  we  have  a  green  in  direction  k 
which  is  interrupted  by  the  green  for  the  extra  phase.  This  can  be  written  as 


(7.2) 


Furthermore,  we  have  that  the  start  of  green  on  the  extra  phase  occurs  after 
the  green  duration  in  direction  k  or 


♦l  +  Sj  +  ’j  *  *it  (7>3) 

This  is,  in  fact,  a  loop  constraint  with  zero  as  the  integer.  Although  we  have 
illustrated  particular  cases,  it  is  in  principle  the  method  for  handling  extra 
phases  or  other  timing  problems  by  insuring  timing,  sequencing  and  by  appropri¬ 
ately  equating  phase  duration. 

We  must  also  specify  input  data  for  an  extra  phase.  Again  we  are 
given  f?  +  f  and  percent  of  through  movement  *  a  and  percent  turning  movements 

such  that  a  +  b  =  1  on  link  1. 
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For  link  i  the  total  flow 'is  fy  *  a(fp  +  f$)  and  for  link  j»fy  *  b(fp  +  f$), 

see(Fig.  7.2.  Link  m  is  handled aspreviously  indicated.  Appropriate  ad justments 
are  made  for  turning  movements »  then  for  trucks  and  buses.  We  shall  have  equal 
platoon  lengths  on  links  i  and  j  and  also  assume  that  link  lengths,  lost  times, 
and  minimum  red  times  are  equal  on  links  i  and  j. 


data 


Since  the  saturation  flow  = 


number  of  lanes 
“  leadway 


we  use  the  following 


Table  7.1 


UTCS-1  Test' Network:  Off  Peak 


Link  Number  Number  of  Lanes  Headway 


117 

2 

2.4 

118 

2 

2.4 

123 

1 

2.4 

124 

2 

2.4 

137 

2 

2.4 

138 

1 

2.4 

141 

2. 

2.4 

142 

2 

2.4 

In  addition  link  113  has  a  25.0  sec.  minimum  red  time. 

If  fs  >  fp  on  link  1  (the  secondary  flow  exceeds  the  primary  flow). 


Figure  7.3  Schematic  for  platoon  lengths  when  f$  >  fp. 


the  platoon  length  on  link  1  then  becomes  (see  Fig.  7.3), 


P  =  9(1  +  fp/ f  j) 

where  g  is  the  green  time  on  link  j.  Further  adjustments  to  this  platoon 
length  are  as  previously  Indicated.  Also  offset  across  link  1  Is  measured  from 
the  beginning  of  the  green  on  j  rather  than  from  the  beginning  of  green  on  link 
k  as  usual. 
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Tables  7.2  thru  7.6  represent  Input  data  and  computational  results  for 
the  UTCS-1  Test  Network:  Off  Peak. 

Table  7.2 


Loops  of  the  UTCS-1  Network  and  Their  Corresponding  Integer  Variables 


Loop  No.  Integer  Variable  _ Link  Numbers 


1 

N1 

101 

■no 

2 

N2 

102 

111 

3 

N3 

103 

112 

. 

4 

N4 

118 

119 

5 

‘  N5 

128 

129 

6 

N6 

'  140  ' 

141 

7 

N7 

130 

131 

8 

N8 

137 

136 

9 

N9 

145 

146 

10 

N10 

113 

114 

11 

Nil 

124 

125 

12 

NI2 

139 

132 

136 

13 

N13 

115 

121 

116 

105 

14 

NI4 

104 

113 

115 

120 

15 

N15 

122 

119 

106 

116 

16 

N16 

125 

120 

126 

131 

17 

N17 

121 

132 

127 

126 

18 

N18 

122 

128 

133 

127 

19 

N19 

141 

133 

145 

139 

20 

N20 

130 

135 

143 

147 

21 

N21 

135 

143 

144 

137 

105 


Table  7.3 


Link  Data  for  UTCS-1  Test  Network  -  Off  Peak* 


Link 

Travel 

Sat. 

Link 

Length 

Time 

Flow 

Flow 

No. 

jftJL 

(sec.) 

(yeh/sec). 

(veh/sec) 

101 

270 

9.82 

.101 

.417 

102 

340 

12.10 

.074 

.833 

103 

430 

16.50 

.112 

.417 

104 

270 

9.82 

.158 

.417 

105 

340 

12.10 

.236 

.833 

106 

430 

16.50 

.203 

.833 

107 

270 

9.82 

.176 

.833 

108 

340 

12.10 

.204 

.833 

109 

430 

16.50 

.221 

.833 

110 

270 

9.82 

.125 

.417 

111  - 

340 

12.10 

.056 

.417 

112 

430 

16.50 

.110 

.417 

113 

220 

10.71 

.163 

.417 

114 

220 

11.70 

.132 

.833 

115 

220 

10.43 

.274 

.833 

116 

220 

11.65 

.253 

.833 

117 

220 

9.30 

.116 

.833 

118 

220 

9.30 

.246 

.833 

119 

220 

10.11 

.265 

1.25 

120 

220 

11.39 

.103 

.833 

121 

340 

15.63 

.157 

.833 

122 

430 

17.33 

.224 

.833 

123 

160 

5.68 

.013 

.417 

124 

160 

5.68 

.013 

.033 

125 

160 

7.33 

.239 

.417 

126 

260 

14.06 

.213 

.833 

127 

260 

12.16 

.368 

.833 

128 

260 

10.18 

.267 

.833 

129 

260 

10.35 

.279 

.833 

130 

200 

10.84 

.251 

1.25 

131 

200 

11.(8 

.338 

1.25 

132 

340 

12.64 

.104 

.417 

133 

430 

18.69 

.242 

.833 

134 

570 

24.57 

.172 

.833 

135 

170 

6.45 

..186 

1.66 

*  See  footnote  on  next  page. 


Platoon 

Length 

J&clel 

.511 

.547 

.502 

.511 

.547 

.525 

.663 

.526 

.579 

.501 

.526 

.579 

.413 

.643 

.474 

.601 

.552 

.552 

.373 

.640 

.723 

.995 

.640 

.640 

.250 

.601 

.626 

.670 

.640 

.447 

.986 

.494 

.424 

.765 

.669 
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Table  7.3  (cont’d.) 


Link 

No. 


136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 
U8 


Link 

Travel 

Length 

Time 

Flow 

Sat.  Flow 

(«.) 

(sec.) 

(veh/secX 

(veh/sec) 

380 

13.51 

.211 

1.25 

380 

13.15 

.173 

.833 

380 

13.15 

.061 

.417 

170 

7.10 

.268 

.833 

370 

11.69 

.289 

.833 

370 

13.51 

.156 

.833 

370 

13.51 

.  .028 

.813 

310 

11.20 

.162 

1.66 

310 

12.36 

.300 

.833 

470 

16.67 

.302 

1.25 

470 

17.53 

.206 

1.2S 

270 

11.70 

.068 

.833 

340 

14.06 

.121 

.833 

Platoon 

Length 

(cycle) 

.375 

.590 

.590 

.542 

.507 

.580 

.580 

.386 

.624 

.682 

.816 

.405 

.480 


Input 

Links 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 
163 


.242 

.075 

.100 

.252 

.279 

.115 

.201 

.455 

.315 

.479 

.345 

.030 

.310 

.134 

.343 


.833 
.833 
.417 
.833 
1.25 
.417 
.833 
i  .25 
1.25 
.833 
.833 
.417 
1.25 
.833 
.833 


*(a)  Lost  time  for  each  link  1  =  4.0  secs. 

(b)  Platoon  length  on  input  links,  p  =  1  cycle. 

(c)  Minimum  red  time  is  15.0  seconds  on  every  link  except  link  113,  where 
it  is  25.0  seconds. 
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Table  7.4 


Results  and  Statistics  for  UTCS-1  Test  Network  -  Off  Peak 


I. 


I 


V 

i  • 


Rows 

674 

Columns 

996 

Variables 

1670 

Integer  Variables 

21 

Non-zero  Elements 

2798 

Density 

.41 

Time 

Iteration 

Functional 

(min) 

No. 

Value 

Continuous  Optimum 

.62 

604 

131.606 

First  Integer  Solution 

1.53 

970 

164.56 

Second  Integer  Solution 

3.68 

1644 

164.56 

Notice  that  the  second  integer  solution  reached  the  same  value  for 
the  objective  function  as  the  first,  tius  the  solution  was  not  improved.  We 
set  a  time  limit  of  four  minutes  of  CPU  time.  Although  the  program  did  not 
prove  optimality  under  this  restriction,  it  seems  that  the  first  Integer 
solution  found,  although  perhaps  not  optimal  is  quite  near  optimal.  Other 
computational  experience  (see  test  network  in  previous  sections)  tends  to 
support  this  hypothesis. 


Table  7.5 


Integer  Nodes  and  Functional  Values 
UTCS-1  Test  Network  -  Off  Peak 


Functional  Value 

164.56 

164.56 

Frequency 
w  =  {ser)“* 

0.01085 

0.01085 

Integer 

Variable 

Integer 

Integer 

N1 

0 

0 

N2 

0 

0 

N3 

0 

0 

N4 

0 

0 

N5 

1 

1 

N6 

0 

0 

N7 

0 

0 

N8 

0 

0 

N9 

0 

0 

N10 

0 

0 

Nil 

Q 

0 

HI  2 

-1 

-1 

N13 

-1 

-1 

N14 

-1 

-2 

N15 

-1 

-1 

N16 

-2 

-2 

N17 

-2 

-2 

HI  8 

0 

-1 

HI  9 

~2 

-2 

N20 

-2 

-2 

N21 

~1 

-1 

Table  7.6 

Results  of  Computations 
UTCS-1  Test  Network  -  Off  Peak 
(C  =  S2  seconds) 


£ 

ft 

Link 

No. 

Offset 

sec (cycle)  • 

Green  Times 

.  sec(cycle) 

Delay 

Delay  Due 
to  SDF 

Degree  of 
Saturation 

10! 

-1.84 

(-.02) 

53.5  (.58) 

3.59 

0 

.418 

ft 

102 

-2.76 

(-.03) 

50.7  (.55) 

5.81 

0 

.162 

§ 

103 

11.05 

(.12) 

48,8  (.53) 

2.6? 

0 

.507 

£ 

104 

-1.84 

(-.02) 

53.5  (.58) 

4.33 

0 

.653 

te 

& 

105 

-2.76 

(-.03) 

50.7  (.55) 

8.29 

0 

.515 

1 

106 

11.05 

(.12) 

48.8  (.53) 

3.59 

0 

.460 

& 

107 

1.84 

(.02) 

50.7  (.55) 

8.66 

0 

.384 

| _ 

108 

2.76 

(.03) 

48.8  (.53) 

6.18 

0 

.46J 

| 

109 

-11 .05 

(-.12) 

53.5  (.58) 

12.44 

0 

.457 

S- 

i 

110 

1.84 

(.02) 

50.7  (.55) 

2.40 

0 

.545 

if 

ft 

111 

2.76 

(.03) 

48.8  (.53) 

5.16 

0 

.253 

i 

112 

-11.05 

(-.12) 

53.5  (.58) 

12.44 

0 

.328 

i 

113 

-18.43 

(-.20) 

67.3  (.73) 

.184 

0 

.535 

114 

18.43 

.20) 

28.6  (.31) 

2.37 

0 

.511 

i 

115 

13.82 

(.15) 

48.8  (.53) 

1.57 

0 

.621 

1 

116 

33.13 

(.36) 

34.1  (.37) 

10.88 

1.172 

.821 

t 

11/ 

23.04 

(.25) 

16.6  (.18)  25.99 

1.209 

.774 

1 

118 

-15.67 

(-.17) 

55.3  (.60) 

11.71 

0 

.492 

119 

15.67 

(.17) 

28.6  (.31) 

2.76 

0 

.684 

i 

120 

-7.37 

(-.08) 

35.0  .38] 

20.46 

0 

.325 

& 

121 

57.14 

(.62) 

24.9  (.27] 

20.83 

0 

.698 

1 

122 

-20.28 

(-.22) 

28.6  (.31] 

31.61 

1.556 

.867 

£ 

123 

71.89 

(.78) 

8.3  (.09] 

38.34 

0 

.346 

& 

124 

44.24 

(.48) 

35.9  (.39] 

19.45 

0 

.040 

125 

7.37 

(.08) 

67.3  (.73] 

0.0 

.769 

.785 

1 

126 

43.40 

(.46) 

35.9  (.39] 

14.29 

1.704 

.656 

i 

127 

11.06 

(.12) 

59.0  (.64] 

0.0 

C 

.690 

1 

128 

48.85 

(.53) 

33.2  (.36] 

22.95 

1.685 

.890 

s 

129 

43.32 

(.47) 

34.1  .37] 

18.53 

1.663 

.905 

I 

130 

21.20 

(.23) 

30.4  (.33 

5.16 

0 

.608 

ft 

131 

-21.20 

(-.23) 

47.9  (.52 

21.29 

0 

.520 

£ 

132 

-16.59 

(-.18) 

57.1  .62] 

13.09 

0 

.402 

| 

133 

20.28 

(.22) 

36.9  (.40] 

.92 

1.366 

.725 

| 

134 

47.93 

(.52 

23.0  .25] 

25.25 

0 

.826 

1 

135 

0.0 

(0.0) 

61.8  (.67] 

.92 

0 

.167 

no 


Table  7.6  (cont’d) 


Link 

Offset 

Green  Times 

No. 

sec(cycle) 

sec(cycle) 

136 

7.37  (.08) 

38.7  (.42) 

137 

-7.37  (-.08) 

58.1  (.63) 

138. 

28.57  (.31) 

22.1  (.24) 

139 

7.37  (.08) 

47.0  (.51) 

140 

14.75  (.16) 

44.2  (.48) 

141 

-14.75  (-.16) 

45.2  (.49) 

142 

23.04  (.25) 

7.4  (.08) 

143 

6.45  (.07) 

39.6  (.43) 

144 

23.96  (.26) 

46.1  (.50) 

145 

-.92  (-.01) 

40.6  (.44) 

146 

.92  (.01) 

.38.7  (.42) 

147 

32.26  (.35) 

16.6  (.18) 

148' 

13.82  (.15) 

44.2  (.48) 

Delay  Due 

Degree  of 

Delay 

to  SDF 

Saturatioi 

3.59 

0 

.402 

6.36 

0 

.330 

21.57 

0 

.610 

3.59 

0 

.631 

1.29 

0 

.723 

18.34 

0 

.382 

33.18 

0 

.431 

0.0 

0 

.227 

5.90 

0 

.720 

18.16 

0 

.549 

20.09 

0 

.389 

10.23 

0 

.454 

0.0 

•  0 

.303 

Input 

Links 


149 

44.2 

(.48) 

19.26 

0 

.605 

150 

12.9 

(.14) 

36.31 

0 

.643 

151 

28.6 

•  31) 

28.94 

1.453 

.774 

152 

31.3 

•34) 

28.02 

1.667 

.890 

153 

28.6 

(.31) 

28.29 

0 

.720 

154 

53.5 

.58) 

12.63 

0 

.475 

155 

53.5 

(.58) 

12.07 

0 

.416 

156 

39.6 

(.43) 

23.87 

•  1.146 

.827 

157 

40.6 

(.44) 

21.20 

0 

.573 

158 

71.0 

(.  77 

6.36 

.084 

.747 

159 

71.0 

.77) 

3.13 

0 

.539 

160 

22.1 

.24) 

29.68 

0 

.300 

161 

47.9 

|.52) 

15.48 

0 

.477 

162 

53.5 

;.58 

10.23 

0 

.277 

163 

53.5 

'.58) 

14.93 

0 

.710 

7.2  Waltham  Arterial 


An  example  of  an  arterial  is  given  which  contains  11  intersections 
on  Main  Street  in  Waltham,  Massachusetts,  Fig.  7.4.  Tables  7.7  and  7.8  repre¬ 
sent  input  data  while  Tables  7.9,  7.10,  and  7.11  give  results.  Space-time 
diagrams  of  the  results  are  given  in  Fig.  7.5.  Part  (a)  of  the  figure  repre¬ 
sents  westbound  flows,  part  (b)  eastbound  flows,  the  difference  being  caused 
by  multiple  phasing  at  some  intersections.  Only  the  pertinent  phases  in  a 
given  direction  are  considered.  Bandwith,  a  geometrical  quantity,  is  not  a 
criterion  in  our  optimization  procedure.  The  settings  provided  by  MITROP 
minimize  total  delay  by  taking  explicit  account  of  all  traffic  movements 
affected  by  the  arterial  signals. 


Tab'.e  7.' 

Loops  for  the  Arterial  and  their  Corresponding  Integer  Variables 
Loop  No.  Integer  Variable  Link  Numbers 


1 

2 

3 

4 

5 

6 
7 
8. 
9 

10 


N1 

104 

107 

N2 

108 

109 

N3 

110 

111 

N4 

112 

114 

N5 

115 

116 

N6 

117 

118 

N7 

119 

120 

N8 

121 

122 

N9 

123 

125 

N10 

101 

103 

Table  7.8 


Link  Data  for  the  Arterial* 

Platoon 


Link  No. 

Link  Length 
(ft) 

Travel  Time 
{sec  ) 

Flow 

(veh/sec) 

Sat.  Flow 
(veh/sec) 

Length 

(cycle) 

101 

810 

23.82 

.106 

.534 

.499 

102 

810 

23.82 

.204 

1.111 

.499 

103 

810 

23.82 

.270 

1.467 

.612 

104 

1105 

32.50 

.182 

.979 

.713 

105 

1105 

32.50 

.052 

.445 

.713 

106 

1105 

32.50 

.047 

.445 

.543 

107 

1105 

32.50 

.175 

.978 

.543 

108 

755 

22.20 

.223 

.890 

.840 

109 

755 

22.20 

1222 

.890 

.798 

110 

800 

23.50 

.201 

.890 

.678 

111 

800 

23.50 

.228 

1.328 

.848 

112 

360 

10.60 

.152 

.801 

,388 

113 

360 

10,60 

.089 

.534 

.383 

114 

360 

10.60 

.287 

.979 

.600 

115 

230 

6.76 

.162 

1.156 

.388 

116 

230 

6.76 

.336 

1.332 

.663 

117 

830 

24.41 

.233 

1.111 

.855 

118 

830 

24.41 

.233 

1.111 

.755 

119 

410 

12.06 

.208 

1.111 

.588 

120 

410 

12.06 

.230 

1.022 

.750 

121 

340 

10.00 

.298 

.979 

.463 

122 

340 

10.06 

.296 

1.068 

.588 

123 

1050 

30.88 

.224 

1.068 

.915 

124 

1050 

30.88 

.083 

.498 

.564 

125 

1050 

30.88 

.193 

.979 

.564 

Input  Links 


126 

.325 

1.289 

127 

.039 

.445 

128 

.138 

.623 

129 

.096 

.712 

130 

.037 

.667 

131 

.011 

.534 

132 

.106 

.979 

133 

.072 

.489 

134 

.125 

.489 

•  135 

.189 

1.156 

136 

.030 

.890 

137 

.088 

.623 

'  138 

.176 

.979 

139 

.159 

1.088 

140 

.086 

1.200 

141 

.089 

.489 

142 

.076 

.489 

*(a)  Lost  time  for  each  link,  1  =  4.0  secs. 

(b)  Platoon  length  on  input  links,  p  =  1  cycle, 

(c)  Minimum  red  time  on  each  link  is  15  secs. 

(d)  Velocity  on  each  link  is  34  ft/ sec. 
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Table  7.9 


Results  and  Statistics  for  the  Arterial 


Rows 

427 

Col umns 

206 

Variables 

631 

Integer  Variables 

10 

Non-zero  Elements 

1725 

Density 

Q.S4 

Time 

Iteration 

Functional 

(min. ) 

No. 

Value 

Continuous  Optimum 

.25. 

410 

79.6803 

First  Integer  Solution 

.53 

571 

91 .8803 

Optimal  Integer  Solution 

.53 

571 

91 .8803 

Optimality  Proved 

.90 

743 

Time  of  Search 

.90 

743 

Number  of  Integer  Variables 

not  at  Continuous  Optimum 

10 

Number  of  Integer  Solutions 

Found 

1 

Table  7.10 


Integer  Nodes  and  Functional  Values  for  the  Arterial 


Table  7.11 


Results  and  Computations  for  the  Waltham  Arterial 
C  =  69  seconds 


Link 

c  Offset 

Green  Times 

Delay  due 

Degree  of 

No. 

sec  (cycle) 

sec  (cycle) 

Delay 

to  SDF 

Saturation 

101 

37.26- 

[.54) 

33.12  (.48) 

6.60 

.414 

102 

37.26 

[.54) 

21.39  (.31) 

6.60 

.592 

103 

44.16 

[-64) 

22.08  (.32) 

10.05 

.575 

104 

46.92 

[.68) 

22.08  (.32) 

14.90 

.581 

105 

73.14 

(1.06) 

14.49  (.21) 

20.32 

.347 

106 

56.58 

[.82) 

13.80  (.20) 

11.84 

.528 

107 

22.08 

[.32) 

48.30  ( 

;.7o) 

0.0 

.255 

1G8 

49.68 

[-72) 

30.36  ( 

[.44) 

13.84 

• 

.569 

109 

19.32 

[-28) 

44.16  ( 

:.64) 

6.20 

.390 

no 

26.22 

[.38) 

44.16  ( 

:.64) 

1.40 

.353 

111 

42.78 

.62) 

19.32  ( 

;.28) 

19.89 

.613 

112 

12.42 

(.18) 

19.32  ( 

[.28) 

11.26 

.678 

113 

35.88 

[.52) 

15.87  ( 

[-23) 

17.34 

1.343 

.725 

114 

-12.42 

[-.18) 

30.36  ( 

:.44) 

8.33 

.666 

115 

19.32 

[.28) 

14.49  ( 

:.2i) 

6.24 

.667 

116 

-19.32 

[-.28) 

48.99  ( 

:.7i 

4.68 

.355 

117 

46.23 

[.67) 

34.50  { 

•  50 

12.39 

.401 

118 

22.77 

[.33) 

53.82  < 

1-78) 

0.0 

.269 

119 

-  1.38 

[-.02) 

53.82  ( 

;.78) 

0.0 

.240 

120 

1.38 

[-02) 

32.43  ( 

:.47) 

12.19 

.479 

121 

9.66 

-14) 

32.43  ( 

:.47) 

0.0 

.648 

122 

-  9.66 

[-.14) 

29.67  ( 

;.43) 

18.68 

.645 

123 

48.30 

[-70) 

29.67  ( 

[-43) 

16.72 

.488 

124 

50.37 

[.73) 

19.32  ( 

[-28) 

9.86 

.595 

125 

20.70 

[,302 

48.99  ( 

:,7i) 

Q.Q  . 

.278  _ 

Inputs 

Links 

126 

18.63  (.37) 

18.77 

.681 

127 

12.42  i 

(.18) 

25.60 

,487 

128 

17.25 

(.25 

•24.10 

2.423 

.886 

129 

14.49 

(.21) 

i  24.36 

.642 

130 

14.49 

(.21) 

i  23.35 

.264 

131 

6.9 

(.10) 

l  28.13 

.206 

132 

11.73 

(.17) 

1  26.19 

.637 

133 

17.25 

(.25) 

1  23.06 

.589 

134 

25.53 

(.37) 

1  18.90 

.480 

.691 

135 

16.56 

(.24) 

1  23.79 

.681 

136 

17.25 

(.25) 

1  21.56 

.135 

137 

13.11 

(.19) 

I  25.41 

.574 

138 

20.01 

(.29) 

)  21.51 

.620 

139 

53.13 

(.77) 

>  •  2.40 

.193 

140 

8.28 

.12)  28.14 

.597 

141 

27.6 

(.40)  16.40 

.455 

142 

35.19 

(.51)  10.93 

.305 

• 
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:  WALTHAM  ARTERIAL  (Main  Street) 
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[:  Time-Space  Diagram  for  Waltham  Arterial,  Eastbound. 


8.  DISCUSSION  OF  RESULTS 


The  preceding  chapters  described  the  model  formulation  of  the  traffic 
signal  network  timing  problem  and  the  computational  aspects  of  applying  mixed- 
integer  linear  prograir.il rig  techniques  (MILP)  for  its  optimization.  The  tech¬ 
niques  were  used  to  calculate  offsets  (Chapter  4)  offsets  and  splits  (Chapter 
5),  and  offsets,  splits  and  cycle  time  simultaneously  (Chapter  6). 

The  conventional  procedure  for  determining  the  signal  control  variables 
is  to  use  a  sequential  decision  process:  A  common  cycle  time  is  selected  for 
the  network  first.  Then  the  splits  at  each  intersection  are  determined  ac¬ 
cording  to  the  proportions  of  demand/capacity  ratios  on  conflicting  approaches. 
Finally,  linking  of  the  signals  is  achieved  by  an  appropriate  method  for  selec¬ 
tion  of  a  fundamental  set  of  offsets  throughout  the  network. 

Experience  of  researchers  and  practitioners  in  the  urban  traffic  control 
field  has  shown  that  cycle  time  may  well  be  the  most  important  control  variable 
for  determining  signal  settings  in  a  network  [38].  Current  approaches  for  se¬ 
lecting  a  cycle  time  can  be  divided  into  two  classes.  The  first  class  is  the 
node  approach.  Since  through  capacity  increases  with  cycle  length,  this  ap¬ 
proach  is  based  on  analyzing  the  capacity  requirements  of  each  intersection  in 
the  network.  The  common  cycle  time  is  determined  according  to  the  requirements 
of  the  most  heavily  loaded  intersection,  i.e.,  the  intersection  with  the  highest 
sum  of  demand/capacity  ratios  on  conflicting  signal  phases.  A  procedure  that 
is  used  for  a  single  intersection,  such  as  Webster's  method  [40],  is  then  used 
to  calculate  the  cycle  length.  This  approach  has  been  primarily  used  in  con¬ 
junction  with  offset  optimization  methods  such  as  Combination  and  TRANSYT  [17]. 
The  main  deficiency  in  this  approach  is  that  the  interaction  of  flows  in  the 
spatial  road  network  structure  of  the  area  is  disregarded.  A  formula  devised 
for  an  isolated  intersection,  assuming  randomly  distributed  arrival  times  of 
cars,  is  not  necessarily  valid  in  a  network  situation  where  flows  are  fed  from 
adjacent  intersections.  The  result  is  generally  a  cycle  time  that  is  too  long, 
thus  causing  excessive  delays. 

The  second  class  is  the  network  approach.  In  this  case  an  attempt  is  made 
to  select  a  cycle  time  that,  while  satisfying  the  capacity  requirements  at  each 
intersection,  is  also  congruent  with  the  particular  network  structure  at  hand. 
Simple  examples  in  this  category  are  the  arterial  progression  schemes  in  which 
a  cycle  that  produces  maximal  bandwidths  is  selected  according  to  distance  and 
speed  data'(e.g.,  [3,24,25]).  The  underlying  principle  is  that  the  optimal  pro¬ 
gression  for  given  block-length  pattern  and  desired  speeds  is  strongly  depen¬ 
dent  on  cycle  time.  In  a  general  network  this  approach  is  principally  used  by 
SIGOP  [37].  A  predetermined  number  of  cycle  times  are  scanned  in  this  method. 
For  each  cycle,  offsets  are  optimized  by  the  OPTMIZ  subroutine  and  performance 
is  evaluated  by  a  coarse  simulation  of  traffic  flow  through  the  network  (the 
VALUAT  subroutine).  The  optimal  set  of  cycle  and  offsets  is  selected  according 
to  the  results  obtained  by  VALUAT.  TRANSYT  indicates  also  the  possibility  to 
iterate  on  cycle  time  in  conjunction  with  its  hill-climbing  procedure  for  off¬ 
set  selection  [34].  However,  computational  considerations  seem  to  rule  out 
this  possibility  in  practice.  Two  deficiencies  of  the  network  approach  in 
SIGOP  are  apparent:  first,  the  offset  optimization  procedure  determines  a  lo¬ 
cal  optimum  rather  than  global  optimum  and,  second,  stochastic  effects  on  link 
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performance  are  ignored.  These  effects  do  not  affect  the  selection  of  offsets 
at  a  fixed  cycle  time,  but  are  of  prime  importance  in  evaluating  a  range  of 
cycle  times.  They  become  pronounced  as  a  signalized  intersection  approaches 
its  capacity  and  in  an  optimal  procedure  would  deter  the  cycle  time  from  as¬ 
suming  values  close  to  the  minimum. 

This  study  models  both  the  deterministic  and  stochastic  effects,  via  the 
link  performance  function  (Chapter  3)  and  the  saturation  deterrence  function 
(Section  5.5),  respectively.  For  the  first  time  a  method  is  presented  that 
optimizes  all  the  decision  variables— cycle  time,  splits,  and  offsets— simul¬ 
taneously.  The  results  obtained  with  this  approach  show  promise  of  substantial 
improvements  in  network  performance  when  compared  to  conventional  procedures. 

The  optimal  settings  for  the  test,  network  in  Chapter  6,  when  all  the  decision 
variables  were  considered  simultaneously  (objective  function  =  53.8010),  showed 
an  improvement  of  11.6%  in  network  performance  with  respect  to  the  settings  ob¬ 
tained  in  Chapter  4,  when  the  decision  variables  were  considered  sequentially 
(objt'tive  function  =  60.8741).  The  interplay  between  the  LPF  and  the  SDF  in 
the  objective  function  for  variable  cycle  time  (for  the  test  network  analyzed 
in  the  report)  is  shown  in  Fig.  6.1  and  is  of  fundamental  importance  in  ana¬ 
lyzing  the  performance  of  area  traffic  control  systems.  Furthermore,  the  re¬ 
sults  in  Fig.  6.5  show  that  the  timings,  and  in  particular  the  cycle  time,  are 
sensitive  to  flows  and  substantial  gains  can  be  achieved  by  setting  the  system 
in  accordance  with  the  actual  traffic  conditions.  From  a  computational  stand¬ 
point  the  method  is  not  at  the  present  time  implementable  as  an  online  real¬ 
time  procedure.  Nevertheless,  MITROP  presents  an  innovative,  comprehensive 
and  systematic  approach  to  traffic  signal  settings  in  a  network  and  should  be 
further  developed  in  several  directions. 

First,  on  an  offline  basis,  the  method  should  be  developed  into  an  opera¬ 
tional  tool  for  use  by  practicing  traffic  engineers.  Results  obtained  in  this 
study  demonstrate  that  sizeable  networks  can  be  optimized  in  relatively  short 
times,  on  the  order  of  a  few  minutes.  Hence,  the  method  Is  capable  of  being 
developed  into  an  Interactive  computer  program  which  can  be  put  on  a  time¬ 
sharing  network  and  made  easily  accessible. 

Second,  the  method  should  be  applied  to  several  urban  networks  to  test 
its  calibration  procedures  and  evaluate  its  performance  on  the  street. 

Third,  the  insight  and  knowledge  gained  in  this  research  can  be  used  for 
modeling  the  traffic  process  in  its  transient  behavior  and  developing  the  ana¬ 
lytical  tools  for  its  online  optimization,  which  is  indeed  the  third-generation 
goal.  However,  there  still  seems  to  be  a  long  way  before  urban  traffic  con¬ 
trol  and  management  can  be  turned  into  an  automatic  feedback  control  system. 

The  traffic  analyst  still  plays  an  important  role  in  the  loop  going  from  traf¬ 
fic  data  to  signal  settings.  Therefore,  first-  and  second-generation  strategies 
will  seemingly  continue  to  constitute  the  backbone  of  traffic  control  procedures 
for  a  long  time  to  come. 
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